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ABSTRACT 



In this thesis a profitability model is developed for 
an exploratory well drilling program. The model considers the 
number of wells drilled, the probability of success, the reserve 
size discovered by a successful well, the discounted profit per 
barrel for oil produced and the cost of an unsuccessful well or 
"dry hole". Each of these parameters are analyzed and a prob- 
ability distribution developed for all but the dry hole costs. 
With this model and the respective probability distributions, 
a FORTRAN IV computer program is developed using the Monte 
Carlo technique of stochastic simulation. The program provides 
the best estimate of discounted profit for the venture, the 
risk of the venture as measured by the standard deviation (and 
the coefficient of variation) and the probability distribution 
of the discounted profit - all after any number of prespecified 
model executions. 

After developing a hypothetical set of base parameter 
data and executing the program, a reference set of values for 
the profitability and risk are obtained. Each parameter set 
is then varied holding the remaining parameters at their base 
values. Using this technique the effects of parameter changes 
on observed model results are assessed. 

Several different sensitivity studies are performed in 
this manner using the data thus developed. These studies 
include the effects of parameter changes on the expected profit 
and risk of the venture using a Significance Ratio technique 
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developed in the thesis. In addition, the convergence of the 
best estimate of discounted profit with the number of model 
executions, as well as the relationship between the number of 
wells drilled and the venture risk are analyzed. 

Also studied is the possible effects of periodicities 
in the combined distribution data with an increased number of 
model executions and the effects of changing a parameter prob- 
ability distribution on the best estimate of discounted profit 
for the venture. 
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CHAPTER I 



INTRODUCTION 



1 . 1 General 

A primary consideration of management in making invest- 
ment decisions is the eventual profitability of the proposed 
venture. There are several ways of measuring profitability 
such as Rate of Return, and Present Worth or Discounted Expect- 
ed Profit. The main purpose in calculating this profitability 
criterion is to summarize into a single measure the quantifi- 
able factors affecting the economic desirability of the project 
under question [25] . 

Profitability calculations have generally been based 
on the best estimates of the individual variables involved. 
However, these best estimates may or may not occur. The re- 
sultant single value of profitability based on these best es- 
timates gives no information on the chances of doing much 
better or much worse then the expected value. 

In his article, Risk Analysis in Capital Investment 
[23], David B. Hertz demonstrates that in certain cases the 
conclusions reached on investment decisions may have been much 
different if the analysis provided not only a single profit- 
ability estimate but also the chances of that estimate being 
exceeded or, on the other hand, not being met. Hertz concludes 
that the application of probabilities will often yield entirely 
different and better decisions. Thus, by developing the pro- 
bability distribution of the profitability, management can 
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make quantitative assessments of the risks involved in approving 
a particular investment. The chances of achieving a minimum 
desired return or of losing money can be determined so that 
management knows the size of the risk it is undertaking. 

In order to determine the profitability of the venture 
and the probability distribution of the profitability, a model 
must be established containing all the factors which bear on the 
problem. In some cases the model is of such a nature that it 
can be solved analytically. That is, the model consists of 
known mathematical procedures or algorithms, which when manip- 
ulated, either directly or iteratively, provide a solution. 

Many types of managerial problems, however, are so complex that 
neither direct nor iterative analytical procedures exist for 
their solution, or, if such methods do exist, they represent 
such great difficulties in application that their utilization 
is impractical. In such cases we turn to the use of a technique 
called simulation. Richards and Greenlaw [42] discuss simula- 
tion as follows: 

Simulation involves designing and utilizing a model 
which replicates some aspect of the firm's operations... 
In general, simulation models may be distinguished from 
other operations research models in that they are non- 
analytical in nature, calling for an experimental app- 
roach, in which workable, satisficing solutions rather 
than optimal ones are usually provided. 

In this thesis an oil exploration investment decision model will 
be developed incorporating several different variables, all of 
which bear on the profitability of an oil exploration venture. 
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Each variable will have a different probability distribution 
making the model so complex that analytical solution is not 
practical. We will therefore turn to simulation techniques, 
namely the Monte Carlo technique (described in Chapter II), to 
provide a measure of profitability and a probability distribu- 
tion of that profitability. 

Before discussing the organization of this thesis it 
may be well to clarify one point which seems to cause a great 
deal of confusion, that is, the distinction between uncertainty 
and risk. Contemporary literature, by and large, still maintains 
F.H. Knight's distinction between these two terms [29]. Risk 
is characterized by situation in which the outcome is not known 
with certainty, but where the probabilities of the outcomes are 
either known or can be estimated. Uncertainty, on the other 
hand, refers to situations for which probabilities of the out- 
comes cannot even be predicted in probabilistic terms. 

Concerning this distinction, Alfred Rappaport [39], 

writes : 



While this distinction may be a useful pedagogical 
device, there is serious doubt about its applicability 
to decision makers in an organizational context. Here 
the distinction is blurred by the fact that decision 
makers generally have some feelings about the probabil- 
ities of future events. Admittedly, these feelings 
range from a degree of confidence bordering on certainty 
to ill-defined feelings approaching mysticism. The 
critical point remains that these feelings affect the 
judgments made by the decision maker. 

Based on the reasons expressed by Rappaport the terms "uncert- 
ainty" and "risk" will be used interchangeably throughout this 
thesis . 
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1 . 2 Organization of the Thesis 

This thesis treats the sensitivity studies in profit- 
ability analysis in two parts. The first part, beginning with 
Chapter II, describes the Monte Carlo method of stochastic 
simulation, a technique used throughout the thesis. In Chapter 
III a general profitability model for an exploratory oil 
drilling venture is developed. Chapter IV discribes some 
general statistical "yardsticks" used to measure central ten- 
dency and variability of probability distributions. Beginning 
with Chapter V and continuing through Chapter VIII each varia- 
ble, or parameter, in the profitability model is analyzed and 
a probability distribution best describing the parameter is 
established. In addition, a hypothetical set of values for 
each parameter is established which will be used as a Base for 
the subsequent sensitivity studies. With the profitability 
model and the respective probability distributions established, 
a FORTRAN IV computer program is developed in Chapter IX using 
the Monte Carlo method to provide the best estimate of dis- 
counted profit for the venture, the probability distribution of 
the profit and the risk of the venture. Finally, in Chapter X, 
using the hypothetical set of values for each parameter pre- 
viously developed, the model is executed and a reference, or 
Base, set of values are obtained for the venture profitability 
and risk. 

The second part of the thesis deals with the actual 
sensitivity studies using the model, parameter probability 
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distributions and Base values developed in previous chapters. 

Chapter XI describes the need for sensitivity analysis 
in general terms and sets the framework for the studies which 
follow. In Chapter XII the results of the sensitivity studies 
are reported and analyzed. Conclusions and recommendations 
resulting from these studies and analyses are contained in 
Chapters XIII and XIV, respectively. 



CHAPTER II 



STOCHASTIC SIMULATION - THE MONTE CARLO TECHNIQUE 
2 . 1 History 

Through the ages man has used chance processes. 

Theories of chance processes were known in ancient times. The 
evidence is in the casting of the astragal bones from the ankles 
of sheep and goats in the manner of dice to subdivide land 
holdings. This was called casting lots. However, it was not 
until the 17th century that probability theory became formalized. 
Interestingly, this mathematical development arose out of in- 
quiries into games of chance. 

Another important development occurred in the 1940 's 
during war work on the atomic bomb. Scientists were frustrated 
by being unable to solve some complex nonprobabilistic mathe- 
matical problems directly by analysis. By combining chance 
processes and probability theory, they solved the problems in- 
directly. To do this, they created a new problem stated in terms 
of stochastic models in such a way that the new problem was equ- 
ivalent to the complex nonprobabilistic problem. But the sto- 
chastic version represented chance processes that could actually 
be carried out. They actually performed the chance processes 
represented by the probabilistic, or stochastic, formulas in 
the new version of their problem. Their methods were equivalent 
to flipping coins or throwing dice, but much more efficient. 

The results of these experiments were then substituted for the 
stochastic expressions of the new version of the problem to give 
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answers to the original nonprobabil istic problem. These scien- 
tists gave the name Monte Carlo to their new technique [5] [22]. 



2 . 2 The Monte Carlo Technique 

Monte Carlo techniques were born when probability theory 
and actual chance processes were used to solve problems that had 
no stochastic aspects at all. Since that time Monte Carlo tech- 
niques have been used to solve a multitude of nonstochastic (or 
deterministic) type problems [33] . 

When dealing with stochastic (or probabilistic) simula- 
tion it is necessary to have a means for dealing with those 
probabilistic variables, the values which assume a frequency 
distribution. Richard F. Barton [4] defines stochastic simu- 
lation as: 



. . . one in which differing outputs trial to trial can 
be obtained without changing the inputs (ignoring random 
numbers as inputs) . Specifically, this means that 
identical parameters, starting conditions, and input 
time path values produce varying outputs trial to trial 
and run to run. 

One commonly used method for working with stochastic variables 

in simulation models is again the Monte Carlo technique. The 

technique is one which both has application to an almost endless 

variety of stochastic simulation problems and at the same time, 

is relatively simple to comprehend and easy to utilize. S.W. 

Hess and H.A. Quigley [24] define the Monte Carlo technique as: 

... a sampling procedure whereby complicated expressions 
involving one or more probability distributions may be 
evaluated. In essence, it consists of simulation an 
experiment to determine some probabilistic property of 
a population of objects or events by the use of random 
sampling applied to the components of the objects or 
events . 
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A parameter value is obtained by drawing randomly from 
the parameter's probability distribution. In a similar 
manner, selections are made for each parameter value 
from its respective distribution. This set of parameter 
values is then substituted into the model and the first 
sample value of the independent variable is computed. 
Subsequent values of the independent variable are obtain- 
ed by drawing additional sets of parameter values. 

The individual values of the independent variable so 
obtained are in approximation to its true probability 
distribution. The approximate distribution approaches 
the true distribution as the sample size increases. 

The Monte Carlo technique actually involves a two sta- 
ge process. In the first stage, stochastic relationships be- 
tween the variables are determined in order to construct the 
simulation model. These relationships may be determined from 
surveys, company data, industry data, executive experience and 
the like. At any rate, once the relationships are determined, 
probabilities of occurrence of frequency distributions are 
assigned for each variable. The frequency distributions are 
then represented by associating groups of numbers to each value 
of the distribution, the size of the group being governed by 
the relative frequency of the value. As an example, assume 
601 of the items in the distribution have a value A, 30% have 
a value B and 101 have a value C. These three values can be 
represented by numbers between 00 and 99. Thus the items 
having a value A would be represented by the numbers from 00 
through 59, those having values of B and C would be represented 
by numbers from 60 through 89 and 90 through 99, respectively. 

Once the model is developed, frequency distributions 
assigned and numbers representing these frequency distributions 
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determined, we move to stage two which involves running the 
Monte Carlo process to obtain a variable value. We do this by 
first selecting a random number (or numbers) from tables if the 
simulation is to be performed manually or by building a random 
number generator into a computerized simulation program. 

(With some computer systems, such as the CDC 6600, it is poss- 
ible to make use of a random number generating function in- 
corporated in the computer systen itself.) Once the random 
number is selected (or generated) it is matched to the corres- 
ponding number assigned to represent each value of the frequency 
distribution of the variable in question. This in turn estab- 
lishes the value of that variable. This process is repeated 
for each variable in the model, the variables combined and a 
single answer is obtained. Stage two is then repeated as many 
times as is necessary to obtain a statistically reliable sample. 
The various stages of the Monte Carlo process have been schem- 
atically diagrammed by M.D. Richards and P.S. Greenlaw [41] 
and are shown in Figure 2.1. 



DETERMINING STOCHASTIC RELATIONSHIP TO CONSTRUCT SIMULATION MODEL 
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CHAPTER III 



THE PROFITABILITY MODEL 



3 . 1 General 

As noted in Chapter I, the profitability of an explor- 
atory well drilling venture and the probability distribution 
of that profitability will be analyzed using a model which 
incorporates the several different factors which bear on the 
problem, namely: the number of wells drilled in the program, 

the probability of a success, the size of field discovered by 
a successful well, the discounted profit per barrel of oil 
discovered and the cost of a dry hole, or an unsuccessful well. 
The model essentially measures the profitability of the venture 
on the basis of Net Present Worth technique whereby the expect- 
ed cash flows are discounted at a suitable discount rate to 
reflect the time value of money. In the model, and throughout 
this thesis, the following notation will be utilized: 

P = Total discounted profit on any one drilling endeavor. 
E(P) = The best estimate (mean) of discounted profit for 
the total exploratory well drilling venture, 
n = Number of wells drilled in the drilling venture, 
p = Probability of a successful well. 

x = Number of successful wells, where x can assume values 
from o to n. 

Q = Reserve size of oil field discovered by a successful 
well (in million bbls) . 
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EQ = Reserve size of oil fields discovered in the venture 
(in million bbls) . 

e = Discounted profit per barrel for oil produced 
(in $ / bb 1 ) . 

C u = Cost for an unsuccessful well, or dry hole. 

Using this notation, the model for a single drilling endeavor 
might be expressed as follows: 

P = x Q e - (n-x) C u (1) 

In an exploratory well drilling venture it must be 
realized that the number of successful wells is a function of 
both the number of wells drilled, n, and the probability of 
success, p. Therefore: 

x = f (n, p) (2) 

Further, the size of reserve discovered in the venture, £ Q, 
is a function of the number of successful wells. 

Z Q = f(x, Q) (3) 

It will be assumed that discounted profit per barrel, e, 
and dry hole cost per unsuccessful well, both as absolute 
values, are independent in that they are not directly derived 
from the number of wells drilled, etc. Therefore: 

P = [zQ = f(x, Q)] e - ^ n - |^x = f(n, p)]j C u (4) 

In Equation 4, which will be used as the profitability 
model in this thesis , it can be seen that the discounted profit 
for the well drilling venture is equal to the size of reserves 
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discovered in the venture (which is a function of the number 
of successful wells and the reserve discovered by each success- 
ful well) multiplied by the discounted profit per barrel for 
oil which will be produced from the reserves, less the number 
or unsuccessful wells (a function of the number of wells 
drilled and the probability of successes) multiplied by the 
cost of the unsuccessful well, or dry hole cost. 

In terms of best estimates, the best estimate of 
discounted profit for the venture, E(P), can be written, 

E(P) = E (x) [ECQ) E(e)] - [n - E(x)] C u (5) 

The relationship expressed in Equation 5, has been used 
by many to estimate the profitability of a venture directly as 
follows : 

Based on all available information, "best guesses" are 
made of the probability of success, which, in this case, is 
the chance that at least one oil field will be found, the 
magnitude of the reserves which would be expected by a success- 
ful well, the profit per barrel for oil produced and the dry 
hole cost. It should be emphasized that normally only one 
value for each parameter is provided, that is the best estimate 
(or guess) . These best estimates are then combined with the 
number of successes determined by 

x = np (6) 

That is, if we had an 18 well program and the probability of 
success was 1 in 9 then we could expect 2 successes. If we 
multiply the 2 successes by the best guess of reserve size per 
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well (say 5 million barrels) we would have a total reserve of 
10 million barrels. If we estimate the most likely discounted 
profit per barrel to be $.80 then we would expect 10 million 
bbls x $.80/bbl or $8 million. From this we would have to 
deduct dry hole costs which would be (18-2) or 16 wells at, 
say, $600,000 per well, or $9.6 million. The best estimate of 
discounted profit for the venture based on this technique is 
therefore - $1.6 million. In all probability the venture would 
not be started on the basis of this analysis. 

As discussed in Chapter I, what this analysis ignores 
is the range of values possible for each of the parameters. 

For example the best estimate for reserve size was 5 million 
barrels per success. That may be the best estimate, but what 
if the reserves actually exceed this best estimate by 251? 

The best estimate of discounted profit would then have been 
$0.8 million and a potentially profitable venture may have 
been passed up. 

In this thesis the basic relationship expressed in 
Equation 4 will be utilized for the profitability analysis. 

In subsequent Chapters each of the parameters will be explored 
to determine their functional relationships and proper frequency 
distributions . 



CHAPTER IV 



DESCRIPTIVE MEASURES OF PROBABILITY DISTRIBUTIONS 



4.1 Measures of Central Tendency (Mean, Median and Mode) 

In analyzing each parameter in the profitability model 
as well as the resultant discounted profit, it is necessary to 
know not only its probability distribution but also some single 
measure typical of the parameter. This single measure will be 
used to describe the central tendency, or the general region of 
continuum in which the distribution lies. 

The most common measure of central tendency is called 
the Mean or Expected Value and is designated x (or x Q ) or E(x). 
The expected value is a measure of location in the sense that 
it roughly locates a middle or "average" value of the variable, 
x [14] . The mean is the sum of a set of values divided by the 
number of values summed. Expressed mathematically the mean of 
a set of n numbers, Xp X 2 , Xp ...., x n is: 

, + x 2 + x 3 + x 4 ... + x n (1) 



E(x) = x = 



n 



Or, where f^ is the frequency of each observation of x (i.e., 
x-p x 2 , Xj, etc.) designated as x^ we can write: 



E(x) = x = 



*i 



n 



( 2 ) 



A second measure of central tendency is the Median . 
The median of a set of values is the middle value when the 
observations are arranged in either ascending or descending 
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order. The median is sometimes defined as the value above and 
below which an equal number of values occur [31] . On a cumu- 
lative probability (or frequency) plot, the median corresponds 
to the value opposite 501 cumulative probability (or frequency) . 



The third measure of central tendency is the mode . The 



mode is defined as the most frequently appearing value. It is 
also referred to as the "most likely value". It should be noted 
that a distribution does not have to have a mode (i.e. the uni- 
form distribution) . Conversely, a distribution may have more 
than one mode (hence the designated "bimodal" for a distribution 
with two modes) . 



of central tendency? Lawrence Mann [32] describes this diff- 
erence as follows: 



This difference is due to the fact that the distribution 
is somewhat skewed. Skewness is the degree to which a 
distribution departs from symmetry. If a distribution 
leans to the right then it is said to be skewed to the 
left (negative skewness); if it leans to the left, then 
it is said to be skewed to the right (positive skewness). 
See Figure 4-1. 



Why is there a difference between these three measures 




x 



X 



Skewed right (positive skew) Skewed left (negative skew) 
Figure 4-1 Illustration of Skewed Distribution 
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In this thesis we will analyze the parameters of the 
profitability model in terms of their mean, median and mode, 
as appropriate. In addition, the expected value and median 
for each case, or parameter set, will be measured. The expected 
value (mean) will be considered the best estimate of discounted 
profit for each parameter set and will be used, along with the 
median, to compare the various cases in the sensitivity studies 
which follow. 



4.2 Measures of Variability (Variance and Standard Deviation) 

The expected value of a random variable is an "average" 
value. However, this measure provides no information about the 
variability of the values of the variable. A measure of this 
variability, or the "spread" or "dispersion" of the values of 
the random variable, adds another dimension to the description 
of the random variable. 

The most useful quantification of uncertainty associ- 
ated with a future event is the standard deviation , o. It 
measures the extent of the spread between individual possible 
outcomes of the event and the expected value of the event. 
Another measure of variability is provided by the variance , 
found by squaring the standard deviation. 

For a discrete random variable, x, the variance is: 



2 = 



(x-iO 2 



n . ? 

z fitx^xr 

_ i 

n- 1 



(3) 
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where x = the mean value of x 
x^ = any value of x 
n = number of values of x 
f^ = frequency of x^ 

The denominator in Equation 2 is replaced by n-1 when the 

statistics are carried out on a sample size n taken from a 

large population - when trying to draw conclusions concerning 

the large population. Division by n-1 in Equation 3 makes 

2 

the sample variance, s x , an "unbiased" estimate of the pop- 

, . . 2 
ulation variance o . 

A 

If x is a continuous variable, the variance is: 

a x 2 = E (x-x) 2 = / (x i - x) 2 g(x) dx (4) 

- 00 



Evaluating the usefulness of variance as a measure of dispersion, 

Samuel Goldberg [15] remarks: 

One difficulty with the variance is that it does not 
measure dispersion in the same units as the values of 
x. Thus, if x has dollar values, then E(x) is a certain 
number of dollars, but since the variance is the mean 
square deviation, Var (x) is measured in dollars 
squared. It is in order to have a measure of disper- 
sion in the same units as the values of x that we 
define the standard deviation as the square root of 
the variance. 

Hence, variance can be adjusted for this disadvantage by taking 
its square root, thus defining the standard deviation, a. This 
adjustment makes the standard deviation directly comparable 
with values of the random variable. The standard deviation 
takes into account, as does the expected value, not only the 
extreme values but all intermediate values, logically relating 
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uncertainty to the expected value [45] . 

In this thesis the standard deviation for the Base and 
all other parameter sets is computed by first calculating the 
variance using a relationship derived from Equation 3 and then 
taking its square root as follows: 



a 



x 





(5) 



4 . 3 Coefficient of Variation 

A knowledge of the expected value and standard deviation 
gives much information about a random variable, even though one 
does not have any knowledge of the individual values of a 
random variable. However, standard deviation is a measure of 
absolute uncertainty, that is, it measures the actual amount 
of variation present in the random variable; and it is depen- 
dent on the scale of measurement. There is a need for additional 
important information to facilitate comparison of variation of 
several random variables [44] . 

The measure which relates uncertainty to the expected 
value of the random variable is the coefficient of variation, 
cv, the formula for which is as follows: 

0 V 

cv = -i (6) 

x 

where cv denotes the coefficient of variation. Use of this 
measure of relative uncertainty makes the variability figures 
directly comparable. It describes the degree of uncertainty 
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or relative risk of the individual expectations about the ex- 
pected value of the random variable [38] . 

Weston and Brigham [47] describe the use of the coeff- 
icient of variation in the following passage: 

... Use the coefficient of variation to compare the 
riskness of alternative investments whenever the stand- 
ard deviation would be misleading. In general, the 
standard deviation is appropriate when percentage rates 
of return are considered, while the coefficient of 
variation is appropriate if dollar returns are consid- 
ered. The reason for this distinction is that rates of 
return are already "standardized" for size of invest- 
ment, but some form of standardization is necessary 
when dollar returns are concerned. 



In the sensitivity studies which follow, the risk 
involved with each Case, or parameter set, will be measured 
in terms of its absolute amount (standard deviation) and in 
terms of its relative value (coefficient of variation) . These 
values will then be related to the standard deviation and 
coefficient of variation of the Base set of parameters indicat- 
ing the change in risk associated with each parameter change. 



CHAPTER V 



PROBABILITY OF SUCCESS 



5 . 1 General 

In assigning the probability of occurrances of future 
events one normally starts by examining past experiences with 
similar events. 

The petroleum industry has been very methodical in 
recording drilling statistics for years. These statistics will 
provide the decision maker with the necessary data to compute 
probability of success, or, as they are called in the industry 
the national success ratios. These national success ratios are 
helpful but may be misleading. They include, for example, all 
wells that produced any amount of oil or gas, regardless of 
whether the amount was small or large, wells drilled in different 
formations, shallow and deep wells, and wells drilled with 
various kinds of skills. 

While the national success ratios provide one set of 
past experiences they may be of little help in assuming prob- 
abilities for a particular venture in a specific location. 

For that particular venture local experiences rather than 
national experience might be more helpful. But even with local 
experience there are still problems in making the statistics 
selective so as to have relevance to the venture at hand. For 
instance, should all wells in all geological provinces be in- 
cluded? Or only those that are producing from the target 
horizon being sought in the venture under consideration? 
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These are difficult questions to answer and may be complicated 
by the fact that if we wanted to exclude all wells not producing 
from a specific geological province a breakdown of data per- 
mitting this exclusion might not be available. 

It should be noted that while local and national ratios 
of successes do provide guides, the assignment of probability 
of success for a particular venture is largely subjective, and 
in some cases, intuitive. But, no matter how the decision is 
arrived at, a probability of success for the venture must be 
stated if a profitability analysis is to be made. Before 
examining how the probability of success will be used in the 
profitability analysis, let us examine some concepts associated 
with probabilities as they relate to business decisions. 

For a long time only interpretation of probabilities 
has been the classical long-run, relative frequency arguement, 
i.e., repetition of the event over and over under identical 
conditions. For the majority of business problems, therefore, 
this concept has been almost useless as a practical guide for 
action. Most business decisions concern singular, or only 
occasional, events that may never be repeated - particularly 
under identical conditions. Therefore, choosing an act because 
it offers the best "long-run average" fits only a limited number 
of business decisions. Yet the "long-run" philosophy persists. 
During his research, Grayson [18] reported the following 
comments by operators: 
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If you drill enough holes, the law of averages will 
take care of you. If you keep playing the averages 
the percentages will hit sooner of later. 

When he questioned them as to what they meant by the words 

"sooner or later", "average", etc. most operators could not 

be specific. 

If we say that the probability of success is 1 in 7 

does this mean that if we drill 7 wells we will have one 

success? The answer is no! It means that if an event, such as 

drilling a well, is repeated over and over, the frequency of 

success will tend to approach the 1 in 7 success ratio. It is 

therefore important to remember that the number of successes 

is dependent on the number of wells drilled as well as the 

probability of success. If we consider the 1 in 7 ratio an 

average success ratio then, as Grayson writes [20] 

... an operator who can drill only a small number of 
wells may expect large variations from the "average". 
The implications for a man with limited funds are clear 
If he had unlimited funds, or even a very large amount 
of funds, he might be able to survive these variations 
early losses or long "runs" of failures - while waiting 
for the "average to come in". But with limited funds, 
an operator may be ruined before the "average" is 
realized. It would be similar in example to the man 
who drowned while crossing a stream with an "average" 
depth of two feet. 
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5 . 2 The Binomial Distribution 

An event that has a certain probability of occurrence, 
p, in a definite number of trials, n, is said to be binomially 
distributed . 

When we drill an exploratory well for oil, one of two 
outcomes is possible. We can have a successful well, which as 
we will see in the next section is defined as a well which 
discovers a field of 1 million barrels or greater, or we can 
have an unsuccessful well, or a'dry hole'. If we let p equal 
the probability of success than 1-p, or q , is the probability 
of a failure. When two exploratory wells are drilled three 
things can happen: both can be successful, one can be success- 
ful and one dry, or both can be dry holes. The chance that 
both are successful is pxp or p , that one is successful and 
1 dry ispxq+pxqor 2pq and that both are dry is q x q 
or q . Using this same logic, if three exploratory wells are 
drilled we can have 3 successes (p 3 ) , 2 successes and one dry 
hole (3p 2 q) , 1 success and 2 dry holes (3pq 2 ) or 3 dry holes 

(q 3 ) . Writing these terms in equation format we have: 



P 3 + 3 P 2 q + 3pq 2 + q 3 = (p+q) 3 (1) 

Thus, if n trials are made each of which results in either 
"success" or "failure", if the probability of success on each 
trial is p and of failure on each trial is q=l-p, and if the 
trials are independent, then the probabilities of all possible 
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results of such a trial are given by the various terms of the 
binomial expansion. 



(p + q) 



n - p n ♦ np n "^q * siaiiV-y 



n! 



x! (n-x) ! 



■p x q n ‘ X 



,n 



( 2 ) 



Where the fourth term in Equation 2 is the general term and 
describes the probability of x successful exploratory wells 
in n independent trials. Hence, the binomial distribution: 



w ( x ) = — tt p x q 

x ! (n-x) ! r M 



,x„n-x 



(3) 



n ^ 

The term x ' i ( n ' -x) T °f ten written as n C x (the number of com- 

binations) . Equation 3 can therefore be expressed as: 



w(x) 



C p x q 
n x F 4 



n-x 



(4) 



Again it is emphasized that Equation 3 provides the 
probability of x successes in n trials only if the probability 
of success, p, is constant from trial to trial and the trials 
are independent. 

To illustrate the use of the binomial distribution in 
predicting successful events let us examine the following 
example . 

Suppose we are planning an exploratory well drilling 
program of 20 wells and we estimate the probability of success, 
p, to be 1 in 7 or say 15%. Then we can also say that the 
probability of failure, q is 1-p, or 851. Using Equation 4, 
we can write the probability of having x success in 20 trials 



as : 
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w(x) = 20 C x (.15) x (.85) 2 °- x (5) 

To determine the probability of 0 successes we would set 



x=0, and. 




w(0) = 


20 C 0 (.15)°(.85) 20 '° 


w(0) = 


^ck)T ( - 15)0( - 85)20 


w (0) = 


( . 85) 20 


w (0) = 


. 0388 or 3.881 



Likewise, the probability of 1 success (x=l) can be written: 



w(l) = 


nw (•15) 1 (.85) 19 


w(l) = 


9 17 (. 15) ( . 85) 19 

1.2165 x 10 1 ' 


w(l) = 


20 C- 15) (.0456) 


w(l) = 


.1368 or 13.68 



A similar analysis could be performed for each success until a 
total of 20 were calculated. The terms could however, be read 
directly from tables [43] . For this example we would have the 
following probability for each possible outcome: 
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TABLE 5-1 

PROBABILITY OF SUCCESS - BINOMIAL DISTRIBUTION 



Number of 

Successes (x) w(x) in percent 



0 


3.88 


1 


13.68 


2 


22.93 


3 


24.28 


4 


18.21 


5 


10.28 


6 


4.54 


7 


1.60 


8 


0.46 


9 


0.11 


10 


0.03 


11 


0.00 

# 


20 


• 

0.00 



5 . 3 Simulation of the Binomial Distribution 

In the Monte Carlo studies to follow we will want to 
simulate the binomial distribution in order to determine the 
number of successful wells for each drilling venture. For 
values of p near 50% and for large n it can be shown that the 
Binomial distribution approaches the Normal distribution [12] . 
The mean value, x Q , and standard deviation, a, are given by: 

x Q = np (6) 



a = /npq 



(7) 
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If x Q corresponding to a cumulative frequency, W(x) = 50%, 
which is a shorthand writing for W(x x^) = 50%, and x Q + a 
corresponding to a cumulative frequency, W(x) = 16% are plotted 
on normal probability paper and a straight line drawn between 
them, the probabilities associated with all other values of x 
could be read directly from the graph. Alternatively, the 
Normal distribution function 



w(x) 



1 




(x-x 0 ) 





( 8 ) 



could be solved analytically to obtain the value of w(x) for 
each success, x. 

If the value of p varies substantially from 50% and if 
the value of n is relatively small can the normal distribution 
be used to simulate the binomial distribution? To answer this 
question we will continue the analysis of the 20 well drilling 
venture that had a probability of success of 15%. From Table 
5-1 a cumulative probability, W(x), can be calculated such that 
W(x <, x^) is the chance that the number of successes will be 
x^ or smaller. 
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TABLE 5-2 

CUMULATIVE PROBABILITY OF SUCCESS - BINOMIAL DISTRIBUTION 



Number of Successes (x^) 


w(x^) in Percent 


(W (x 


£ x^) in Percent 


0 




3.88 




3.88 


1 




13.68 




17.56 


2 




22.93 




40.49 


3 




24.28 




64.77 


4 




18.21 




92.98 


5 




10.28 




93.26 


6 




4.54 




97.80 


7 




1.60 




99.40 


8 




0.46 




99.86 


9 




0.11 




99.97 


10 




0.03 




100.00 


11 

• 




0.00 




100.00 


• 

20 




0.00 




100.00 


The values of (Wx 


£ x^) corresponding 


to values of x 


on Table 5-2 were 


plotted 


on normal probability paper (Figure 


5-1) as circular 


points . 


If these values, in 


fact 


approached 


a normal distribution then 


these points would 


fall 


along a 



straight line passing through x Q at (Wx £ x^) = 501 and x q + o 
at W (x < x.) = 84%. 

l 



From Equation (6) 
x Q = np 
x 0 = 2 0 ( . 1 5) 
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From Equation (7) 



o = / npq 

a = / ( 20 ) (. 15 ) (. 85 ) 
a = / 2.55 
a = 1.60 

Therefore, x Q + a = 4.60 



A straight line is drawn through these two points on Figure 5-1. 

It will be noted from Figure 5-1 that the binomial 
distribution points lie off the normal distribution line. It 
can therefore be seen that for small values of p (in the range 
of 15%) and for relatively small values of n, the use of the 
normal distribution would not be suitable to simulate the bi- 
nomial distribution for the entire range of x^. 

With this in mind the binomial distribution for this 
study can be simulated indirectly by assigning, for each value 
of x, a series of consecutive numbers equal to the values of 
w(x) listed on Table 5-1. 



TABLE 5-3 



SIMULATION OF BINOMIAL DISTRIBUTION - RANDOM NUMBERS 



Number 



of Successes (x) 


w(x) x 1000 


Random ] 


Number 


0 


0388 


0000 


- 0387 


1 


1368 


0388 


- 1755 


2 


2293 


1756 


- 4048 


3 


2428 


4049 


- 6476 


4 


1821 


6477 


- 8297 


5 


1028 


8298 


- 9325 


6 


0454 


9326 


- 9779 


7 


0160 


9780 


- 9939 


8 


0046 


9940 


- 9985 


9 


0011 


9986 


- 9996 


10 


0003 


9997 


- 9999 
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If we now generate a four digit random number from some random 
number generator this number can be equated to a specific value 
of x. See Table 5-3. As an example, if a random number of 
9395 were generated this would be equivalent to 6 successes. 
Likewise, if a number 0035 were generated this would equate to 
0 successes, or all dry holes. With enough trials the binomial 
distribution can be simulated in this manner. 



5 . 4 Coefficient of Variation 

From Equations 6 and 7 the coefficient of variation can 
be derived as follows: 



CV 




/ . n .PiL 
np 



(9) 



Since q = (1-p), Equation 9 can be written: 



CV 



/ np (1-p) 
np 



From which we obtain: 




( 10 ) 



Using the data for number of wells drilled and prob- 
ability of success previously noted, it is possible to calculate 
the Coefficient of Variation using Equation 10 as follows: 



CV = 0.53 



(ID 
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CHAPTER VI 



OIL RESERVOIR SIZE 

It goes without saying that the economic success or 
failure of a drilling venture may well be determined by the 
size of the field discovered as a result of a successful well. 
Of course the size of the field is unknown prior to drilling. 

It is therefore essential, if we are going to attempt to es- 
tablish the relative profitability of the venture, to determine 
two items. First, whether there is a general mathematical ex- 
pression which can be used to characterize the probability 
distribution of field sizes in the geological province. In 
other words, do the field sizes vary in such a way that their 
probability of occurrence can be predicted? Secondly, if a 
general field size distribution function can be established as 
a model, what are the boundary conditions we should establish 
for our particular venture? Each of these questions will now 
be explored. 

6 . 1 Field Size Distribution 

At first glance it would seem quite unlikely that a 
general mathematical model of field size distributions could be 
established - considering the complex process of oil reservoir 
formation, or, if by chance a model could be established for 
one geological province it would have general application in 
another province. 

Some of the earliest work in this area was reported by 
J.J. Arps and T.G. Roberts [3] in 1958 after a study of reserve 
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data from the Denver- Julesburg Basin. In their study, a 5.7 
million acre sample was chosen on the east flank of the Basin 
(located in Nebraska, Wyoming and Colorado) and the reserve 
estimates of all 338 fields found to that date in the Lower 
Cretaceous sands were investigated. Arps and Roberts arranged 
the fields in nine groups according to their areal extent and 
estimated the ultimate primary recovery from each group. The 
limits of these groups were chosen in such a way that each 
group represented fields about twice as large as the preceeding 
one. The nine groups thus arranged are represented on Table 6-1. 

TABLE 6-1 

OIL RECOVERY DATA ON DENVER- JULESBURG BASIN FIELDS 





Range of Productive 


Number of 


Estimated Average 


Group 


Acres 


Fields 


Ultimate Recovery 
(M bbls) 


1 


28 


56 


10.1 


2 


28 - 57 


77 


53.0 


3 


57 - 113 


41 


165.0 


4 


113 - 226 


70 


295.7 


5 


226 - 453 


40 


957.6 


6 


453 - 905 


33 


1,876.2 


7 


905 - 1810 


19 


4,481.6 


8 


1810 - 3620 


1 


11,460.0 


9 


3620 


1 


50,750.0 



When a frequency - density distribution of the data 
presented in Table 6-1 was plotted versus Average Ultimate 
Recovery on logarithmic-probability paper a straight line can 
be best fit through the resulting points. (See Figure 6-1). 
This demonstrates that the field sizes in the Denver- Julesburg 



Basin are Lognormally distributed. 
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While a Lognormal distribution best fits the data in 
the Denver- Julesburg Basin there is still no reason to general- 
ize and say that all reserves follow a Lognormal frequency dis- 
tribution. If, however, the Lognormal distribution remains 
invariant under a wide variation in the definition of reserves 
and geological regions then a generalization can be made. 

To test this invariancy Gordon M. Kaufman [27], in an 
extensive study reported in 1962, accumulated data from several 
different sources including: The Oil and Gas Conservation 

Board of Alberta Report to the Lieutenant-Governor- in-Council 
as summarized in The Financial Post; Oil and Gas Journal stat- 
istics on two states for 1946 and for 1960; an individual 
company's estimate of reserves in four basins; and the Arps - 
Roberts article on the Denver- Julesburg Basin. 

The data accumulated by Kaufman, covering fields in 
four coastal basins, fields in North and South Louisiana (both 
in 1946 and in 1959), fields in Oklahoma, etc. were plotted on 
sixteen (16) separate graphs on logarithmic-probability paper. 
Kaufman, in analyzing the results of these plots, writes [28] 

The graphs displaying this data tentatively suggest 

that : 

1. Within a given area the functional form of empirical 
histograms of reported field sizes does not change 

(a) over time, 

(b) with changes in the definition of "reserves", 

(c) with changes in the minimum size of field 
reported, or 

(d) from one geographic area to another. 

2. The* functional form is the same for a well-defined 
geological basin such as the Denver- Julesburg as for 
an arbitary geographic area such as a state. 
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Although each of /the sixteen (16) graphs^7 is a plot 
of data from field's of varying geological ages and of 
varying length of production lives, the invariance in 
the functional forms is striking enough to warrant the 
conjecture that added data blocked on these two factors 
would not alter the conclusions stated above. 



A visual observation of these graphs reveals that most points 
generally fall along a straight line indicating a Lognormal 
frequency distribution of field size. 

In another analysis of oil field size distributions, 
Folkert Brons reports the result of a study of the predicted 
ultimate recovery of 75 Southern Louisiana Miocene Fields [10]. 
As shown on Figure 6-2, when the data is plotted on logarithmic- 
probability paper the points fall along a straight line indi- 
cating a Lognormal distribution function for these fields. 

On the basis of the evidence presented in the studies 
cited above, of reservoirs with varying reserve sizes, geologic 
origin and physical characteristics, it will be assumed for 
further analysis that a Lognormal distribution function can 
serve as a general mathematical model for reservoir size 
distributions . 

6 . 2 The Lognormal Distribution 

The lognormal distribution can be defined as the dis- 
tribution of a random veriable whose distribution obeys a 
normal distribution function. Thus, if x is lognormally dis- 
tributed then y = In x is normally distributed and, 



w(y) 




CD 
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10 100 1,000 
Ultimate Recovery, Million Barrels 
Figure 6-2 Ultimate Recovery - Southern Louisiana Miocene Fields 
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If the distribution is desired in terms of x rather 
than y a change in variables is performed which results in the 
probability density function for y expressed in terms of x as 
follows : 



w(x) 




( 2 ) 



x 



where: x = independent variable 
Cy = standard deviation 
y 0 = mean value of y 

The equation gives a lognormal distribution curve, 
with its characteristic positive skew, as shown below. 




Lognormal Distribution - Positive Skew 



Figure 6-3 
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Since y is normally distributed: 
y m (median = y Q (mean) = y(mode), 



( 3 ) 



but since x = e^, x is not normally distributed and the values 
of x m (median), X Q (mean) and x(mode) must be calculated separ- 
ately [9] : 

The Median, x m : 



y = In x 

y. 



x m = e 
m 



m 



(4) 



since y m = y Q 



x = e 
m 



and, In x m = y Q 



(4A) 

(4B) 



The Mean, x Q : 




X 

o 

II 

CD 

o 


(5) 


i + °y 2 

or, In x Q = y Q + 


(5A) 


The Mode, x mode : 




y 0 2 

x mode e y 


(6) 


or, In x moc j e yo ” a y 


(6A) 
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2 

Further, the Variance , a x is given by the following expression 



a v 2 



a x 2 = exp { 2y 0 + a 2 } (e ^ -1) 



(7) 



Inserting Equation 5 in Equation 7 we obtain: 

a x 2 = x o 2 (e^ 2 " 1 ) (8) 

From which the Coefficient of Variation , CV, may be expressed 
as : 

(9) 



CV = — = e°y 2 -l 
x o 

For small values of Oy we can approximate Equation 9 as: 



CV = a. 



( 10 ) 



6 . 3 Boundry Conditions 

If we assume that the field (or reservoir) sizes in 
any area are lognormally distributed and will fall along a 
straight line when plotted on logarithmic-probability paper 
then the next problem is one of specifing the position of the 
line which would represent the reservoir size distribution in 
the specific location where the drilling venture is to be 
conducted. 

This line could be drawn if someone could predict the 
median value of the reserve size to be discovered by a (yet 
undrilled) successful well (i.e. the reserve size corresponding 
to a 501 probability) and the standard deviation of the dis- 
tribution representing the logarithm of the field sizes. 
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Then, by inserting these values in Equations (4B) , (5A) and 
(6A) previously noted for median, mean and mode a straight line 
could be drawn on log-probability paper. The problem is, how- 
ever, that the operator (or his agent) would probably not be 
able to specify either the median or the standard deviation. 

What he most likely would be able to specify is his estimate of 
the highest and lowest values of reserves he would expect with, 
say, 90% certainty. He would probably base these estimates on 
either his own experience or the experience of others who have 
drilled in the area. As an example, the operator may not be 
able to specify an exact reserve size, although he might say 
that he would expect the reserves to run from 5 million to 150 
million barrels with 90% confidence. If these two points were 
plotted on log-probability paper - with the 5 million barrel 
figure plotted at the 5% probability point and the 150 million 
barrel figure plotted at the 95% probability point a straight 
line could be drawn between them which would represent the 
Lognormal distribution of all reserves in the area. 

With the line plotted, the median value of Q (or x m in 
Equation 4B) , corresponding to a 50% probability, and the value 
of Q median plus one standard deviation, corresponding to an 
84% probability, can be read directly. In addition, to use this 
model in our determination of the relative profitability of the 
venture using the Monte Carlo technique, it is essential that 
we be able to determine a reservoir size and corresponding value 
of probability (that the reservoir size will not be exceeded) 
for each point along the line. 
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To illustrate this procedure, as well as to generate 
data for use in the Monte Carlo studies conducted as part of 
this thesis, let us assume we are going to drill for oil on- 
shore in the Southern Louisiana area. By going to the approp- 
riate issue of the Oil and Gas Journal [48] we would find a 
listing of the sizes of all fields in the Southern Louisiana 
on-shore area. These fields could then be arranged in ascend- 
ing order of size (Q) and a relative frequency density, w, and 
cumulative frequency, W(Q < Q^) calculated. This tabular 
listing, along with the calculation of relative and cumulative 
frequency is shown as ACTUAL Values in Table 6-3. 

Having the actual values of cumulative frequency cal- 
culated, these values can be plotted on log-probability paper 
with the ordinate being the cumulative frequency and the ab- 
scissa being the field size. The triangular points on Figure 
6-4 represent a plot of the actual data for the Southern Louis- 
iana on-shore area. 

The operator can now use these points (or any other 
data he has available) as a guide in drawing a straight line 
on the log-probability paper beginning at some upper value of 
Q, corresponding to a 951 probability, and at some lower value 
of Q, corresponding to a 5% probability. In this example the 
points chosen were 79 million barrels as the upper value and 
5 million barrels as the lower value. With this line drawn 
(Curve A, Figure 6-4) it can now be said that, for the purposes 
of this profitability analysis, 901 of all reserves discovered 
by a successful well will fall along the straight line, 51 of 
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them will exceed the upper value and 5% will be less than the 
lower value. Theoretically then, there is a 5% possibility 
that a reservoir could be discovered with infinite reserves on 
one extreme and a 5% possibility of a reservoir with a 0 res- 
erve on the other. Practically of course, this is not reason- 
able since 0 reserve means no oil and an unsuccessful well. 

If this straight line represents the reserve size distribution 
for a successful well then some lower limit must be established 
which would represent the reserve size below which the well in 
question would not be considered successful, even though some 
oil was discovered. In this example a lower limit of 1 million 
barrels was chosen. Therefore, we say to be considered succ- 
essful a well must discover at least a 1 million barrel reserve. 

At the upper extremes it would be foolish to say that 
an infinite reservoir could be discovered even though this is 
possible by an extrapolation of Curve A, Figure 6-4, to a 1001 
probability point. Realizing that this is not physically 
possible an arbitrary upper limit should be established. In 
this example it was observed that in the Southern Louisiana 
on-shore area the maximum field size is 137 million barrels 
(see Table 6-3, ACTUAL). On this basis an upper limit of 140 
million barrels was chosen. 

6 . 4 The Lognormal Distribution of Field Sizes 

With the boundary conditions now established and Curve 
A, Figure 6-4 drawn, it is now possible to proceed with the 
calculation of the mean, median and mode as follows: 
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The Median, x m : 



Q at 50% = 20.0 million barrels 



y Q = In 20.00 

y = 2.996 
1 o 

Y 0 +a y = In x (at W ^ x = 84%) 
y 0 + ° y = In 46.0 
2.996 + o y = 3.829 
= 0. 833 




From Equation 4B 




m 



Oy. 2 = 0.694 



a y 2 = 0.347 



2 



The Mean, x Q : 



From Equation 5A 



? 




In x 



2.996 + 0.347 



o 



In x 



o 



3.343 



x 



28.3 



million barrels 



o 
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The Mode, x mode : 

From Equation 6A 

ln x mode y o a y 2 
In x mo( j e = 2.996 - 0.694 

ln x mode ’ 2 ' 302 

In x mo( j e = 10.0 million barrels 



From Equation 9 

cv = ^ = \/e a y 2 -i 

x o 

CV = \je®~‘ 694 - 1 

CV = /2. 0017-1 = /I .0017 

CV = 1.00 

If we assign random numbers to represent W(Q <, Q^) we 
can, from Curve A, Figure 6-4, manually determine the corres- 
ponding value of Q as follows: 
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TABLE 6-2 

FIELD SIZES - RANDOM NUMBERS 

W(Q < Q t ) Q 

Random No . (in \ ) (x 10 6 bbl) 



05 


05 


5.0 


06 


06 


5.5 


07 


07 


5.8 


93 


• 

93 


69^0 


94 


94 


73.0 


95 


95 


80.0 



While the process of manually reading the graph, tab- 
ulating the values of Q corresponding to Random Variables 
which represent W(Q <_ Q^) and inputing these numbers as data 
in the Monte Carlo computer program is possible, the process 
is quite laborious. To equip ourselves for an efficient sim- 
ulation system a Lognormally Distributed Variable Generator 
can be used to generate the values of Q on a high speed 
digital computer. 

Such a Generator was initially developed by McMillan 
and Gonzalez [34], A version of the Generator, modified for 
use in this problem, incorporating the concepts of upper and 
lower limits, is included in Appendix A. Using this Generator 
it is necessary only to input the median Value of Q and the 
value of Q at W(Q £ Q^) = 84%, both easily obtainable from 
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the plot of a straight line connecting the 51 and 95% points, 
and the upper and lower limits, as previously discussed. 

With this data inputed, the Generator will provide field sizes 
(Q) which are lognormally distributed within the boundary 
conditions established for the program. 

In order to test the reliability of the deviates 
provided by the Lognormal Generator, 1000 values of Q were 
generated using the program included in Appendix A and the 
boundary conditions previously discussed. The results of this 
test, along with a calculation of relative and cumulative 
frequency, is included on Table 6-3 under the heading of 
SIMULATION. For comparison purposes the values provided by 
the Generator were placed beside the corresponding actual 
values. Additionally, Curve B, Figure 6-4 was plotted with 
points generated by the Lognormal Generator. 

By comparing Curves A and B and the actual data points 
on Figure 6-4 it will be seen that excellent correlation 
exists except in the lower values of Q - but even there the 
difference is not considered significant. On this basis, the 
Lognormally Distributed Variable Generator will be used to 
simulate lognormally distributed values of field size in this 
thesis . 
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Figure 6-4 

Cumulative Distribution of Field Sizes 
Southern Louisiana On-Shore 
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TABLE 6-3 

SIZE OF FIELD - SOUTHERN LOUISIANA ON-SHORE 
(1970 Estimates) 





ACTUAL 


SIMULATION 


Field 


Freq . 


Rel . Freq . 


Cum. Freq . 


Freq . 


Rel . Freq . 


Cum. Freq . 


Size -Q 




f 






f 




A 


(£) 


(w=ffxl00) 


(W < Q) 


(£) 


(w=irxl00) 


(W < Q) 


(xlO bbl) 




(%) 


(%) 




(%) 


(%) 


1 








3 


0.3 


0.3 


2 








11 


1.1 


1.4 


3 








18 


1.8 


3.2 


4 








25 


2.5 


5.7 


5 








24 


2.4 


8.1 


6 


2 


2.4 


2.4 


33 


3.3 


11.4 


7 


3 


3.5 


5.9 


27 


2.7 


14.1 


8 


2 


2.4 


8.3 


34 


3.4 


17.5 


9 


0 




8.3 


19 


1.9 


19.4 


10 


5 


5.9 


14.2 


45 


4.5 


23.9 


11 


4 


4.7 


18.9 


34 


3.4 


27.3 


12 


5 


5.9 


24.8 


25 


2.5 


29.8 


13 


0 




24.8 


23 


2.3 


32.1 


14 


5 


5.9 


30.7 


37 


3.7 


35.8 


15 


7 


8.2 


38.9 


28 


2.8 


38.6 


16 


3 


3.5 


42.4 


37 


3.7 


42.3 


17 


4 


4.7 


47.1 


32 


3.2 


45.5 


18 


5 


5.9 


53.0 


25 


2.5 


48.0 


19 


2 


2.4 


55.4 


25 


2.5 


50.5 


20 


1 


1.2 


56.6 


27 


2.7 


53.2 


21 


1 


1.2 


57.8 


21 


2.1 


55.3 


22 


1 


1.2 


59.0 


18 


1.8 


57.1 


23 


4 


4.7 


63.7 


22 


2.2 


59.3 


24 


1 


1.2 


64.9 


20 


2.0 


61.3 


25 


0 




64.9 


11 


1.1 


62.4 


26 


1 


1.2 


66.1 


19 


1.9 


64 . 3 


27 


1 


1 . 2 


67 . 3 


15 


1.5 


65.8 


28 


3 


3.5 


70.8 


12 


1.2 


67.0 


29 


1 


1.2 


72.0 


13 


1.3 


68.3 


30 


1 


1.2 


73.2 


20 


2.0 


70.3 


31 


1 


1.2 


74.4 


15 


1.5 


71.8 


32 






74.4 


8 


0.8 


72.6 


33 






74.4 


13 


1.3 


73.9 


34 






74.4 


18 


1.8 


75.7 


35 


2 


2.4 


76.8 


6 


0.6 


76.3 


36 


1 


1.2 


78.0 


13 


1.3 


77.6 


37 






78.0 


4 


0.4 


78.0 


38 






78.0 


8 


Q . 8 


78.8 


39 






78.0 


7 


0.7 


79.5 


40 


2 


2.4 


80.4 


6 


0.6 


80.1 



ie 

iz 

0 6 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 
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TABLE 6-3 Continued 



ACTUAL 


SIMULATION 


Freq . 


ReLFreq . 


Cum. Freq . 


Freq . 


Rel . Freq . 


Cum. Freq . 


(f) 


(w=jrxl00) 

(%) 


(W < Q) 

f%) 


(f) 


£ 

(w==rxl00) 

m 


(W < Q) 
(%) 


1 


1.2 


81 . 6 


5 


0.5 


80.6 


2 


2.4 


84.0 


8 


0.8 


81.4 






84.0 


9 


0.9 


82.3 






84.0 


3 


0.3 


82.6 


1 


1 . 2 


85.2 


6 


0.6 


83.2 






85.2 


8 


0.8 


84.0 


1 


1.2 


86.4 


4 


0.4 


84.4 






86.4 


9 


0.9 


85.3 






86.4 


4 


0.4 


85.7 






86.4 


2 


0.2 


85.9 


1 


1.2 


87.6 


9 


0.9 


86.8 






87.6 


5 


0.5 


87.3 


1 


1.2 


88.8 


5 


0.5 


87.8 






88.8 


2 


0.2 


88 . 0 






88.8 


6 


0.6 


88.6 






88.8 


5 


0.5 


89.1 






88.8 


2 


0.2 


89.3 






88.8 


7 


0.7 


90.0 


1 


1 . 2 


90.0 


3 


0.3 


90.3 


2 


2.4 


92.4 


6 


0.6 


90.9 






92.4 


6 


0.6 


91.5 






92.4 


2 


0.2 


91.7 






92.4 


7 


0.7 


92.4 


1 


1. 2 


93.6 


2 


0.2 


92.6 






93.6 


1 


0.1 


92.7 






93.6 


5 


0.5 


93.2 






93.6 


1 


0.1 


93.3 






93.6 


2 


0.2 


93.5 






93.6 


4 


0.4 


93.9 






93.6 


3 


0.3 


94.2 






93.6 


3 


0.3 


94.5 






93.6 


3 


0.3 


94.8 






93.6 






94.8 






93.6 


3 


0.3 


95.1 


1 


1.2 


94.8 


1 


0.1 


95.2 






94.8 


2 


0.2 


95.4 






94.8 


1 


0.1 


95.5 






94.8 


1 


0.1 


95.6 






94.8 


1 


0.1 


95.7 






94.8 






95.7 






94 . 8 


1 


0.1 


95.8 






94.8 


3 


0.3 


96.1 



Fie 

Siz 

: 10 6 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 



52 



TABLE 6-3 Continued 



ACTUAL 


SIMULATION 


Freq . 


Rel . Freq . 


Cum . Freq . 


Freq . 


Rel . Freq . 


Cum . Freq . 


(£) 


(w=|xl00) 


(W < Q) 


(f) 


(w=£xl 00) 


(W < Q) 




(*) 


(%) 




(%) 


(%) 






94.8 






96.1 






94.8 


1 


0.1 


96.2 






94.8 


4 


0.4 


96.6 






94.8 






96.6 






94.8 






96.6 






94.8 


2 


0.2 


96.8 






94.8 






96.8 






94.8 


1 


0.1 


96.9 






94.8 






96.9 






94.8 






96.9 


l 


1 . 2 


96.0 


1 


0.1 


97.0 






96.0 


3 


0.3 


97.3 


l 


1 . 2 


97.2 


4 


0.4 


97.7 






97.2 


2 


0.2 


97.9 






97.2 






97.9 






97.2 






97.9 






97.2 






97.9 






97.2 






97.9 






97.2 


2 


0.2 


98.1 






97.2 


1 


0.1 


98.2 






97.2 


1 


0.1 


98.3 






97.2 






98.3 






97.2 


1 


0.1 


98.4 






97.2 


1 


0.1 


98.5 






97.2 






98.5 


l 


1 . 2 


98.4 






98.5 






98.4 






98.5 


l 


1.2 


99.6 






98. 5 






99.6 






98.5 






99.6 






98.5 






99.6 


1 


0.1 


98.6 






99.6 


1 


0.1 


98.7 






99.6 






98.7 






99.6 






98.7 






99.6 


1 


0.1 


98.8 






99.6 






98.8 






99.6 


1 


0.1 


98.9 






99.6 






98.9 






99.6 






98.9 






99.6 






98.9 






99.6 


3 


0.3 


99.2 






99.6 






99.2 






99.6 


1 


0.1 


99.3 
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TABLE 6-3 Continued 



Field 
Si ze-Q 

(xl0 6 bbl) 


ACTUAL 


SIMULATION 


Freq. Rel.Freq. Cum.Freq. 

(f) (w=|xl00) (W < Q) 

(%) (%) 


Freq. Rel.Freq. Cum.Freq. 

(f) (w4xl00) (W < Q) 

m m 


126 


99.6 


99.3 


127 


99.6 


99.3 


128 


99.6 


99.3 


129 


99.6 


99.3 


130 


99.6 


99.3 


131 


99.6 


99.3 


132 


99.6 


99.3 


133 


99.6 


99.3 


134 


99.6 


1 0.1 99.4 


135 


99.6 


99.4 


136 


99.6 


99.4 


137 


1 1.2 100.0 


99.4 


138 


100.0 


1 0.1 99.5 


139 


100.0 


99.5 


140 


100.0 


5 0.5 100.0 
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CHAPTER VII 



DISCOUNTED PROFIT PER BARREL 



7 . 1 General 

In a very real sense, the ultimate profitability of an 
exploratory well drilling venture, where oil is discovered, is 
dependent not only on the size of the new reserves but also 
(and perhaps more importantly) on how much profit the operator 
will receive for each barrel of oil produced. The case can be 
visualized where the operator in a successful exploratory drill- 
ing venture discovers vast potential oil reserves but the 
expence of production far exceeds the posted price he can 
receive for the oil produced. In this case no matter how large 
the potential reserve, the profitability of the venture is less 
than 0, or if produced, the operator would suffer a loss. The 
cases of off-shore drilling in deep water and production on the 
North Slope of Alaska (the latter receiving wide publicity) 
illustrate this notion. It is therefore essential in a profit- 
ability analysis to estimate the expected profit per barrel on 
reserves discovered by a successful well. In practice, though, 
it is rare that the operator would commit himself to "the" 
expected profit per barrel for the venture. He could perhaps 
provide his estimates of the lowest, highest and most likely 
values of profit per barrel basing these estimates on his 
experience in the area or on factors similar to those which 
permitted him to estimate the probability of success (see 
Chapter V) . 
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7 . 2 Discounted Profit Per Barrel Estimate 

In estimating the highest lowest and most likely values 
of profit per barrel the operator would first have to provide 
his estimate of the highest, lowest and most likely posted 
price for the oil. From this he would have to deduct his ex- 
penses including, royalties, severance taxes, operating costs, 
drilling and completion costs, etc. This would then provide 
his net profit per barrel for each or the three estimates. 

But, since reserves cannot be produced all at once, the profit 
per barrel can not be figured on the basis of current sales 
prices. Production will be spread into the future (say 10 years) 
and future money is not as valuable as present day money. To 
account for this financial fact-of-life the operator will have 
to convert his net profit per barrel calculated at today's 
prices into a discounted cash flow. In this discounting (or 
deferment) process a discounting factor is developed by calcu- 
lating the Present Worth of a $1 annuity (i.e. constant income 
or expenditure) received each year for a specified length of 
time, n, at a specified discounted rate, i, [16]. This is given 
by the expression: 



P = R 



(1 + i) n -1 
1(1 + i) n 



( 1 ) 



where: P is the Present Worth 

R is the amount of the annuity (taken as 



$1 in this development) 
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1 



is the uniform series Present Worth 
or Present Value Factor. 

The Present Worth Factor can also be written as 
or Values of the Present Worth can be 



(1 + i) n -1 
1(1 + i) n 



(1 + D n 



1 



obtained from Tables. If these values are then divided by the 
number of discounting periods, n, the deferment or discounting 
factor is obtained. 

To illustrate, the discount factor based on a discount 
rate of 10% for a period of 10 years is calculated as follows: 
From Equation 1, 



P 



R 



(1 + i) n -1 
i (1 + i) n 



Let R = $1 
i = 10% 
n = 10 years 

From Tables [17] the Present Worth Factor is 6.1444 

P = 1(6.1444) 

P = 6.1444 



Discount 

Discount 

Discount 



factor 

factor 

factor 



= P 
n 

_ 0.144 
10 

= .6144 



( 2 ) 



Concerning this technique, Grayson [19] writes: 

Many operators, ... object to the present value concept 
First, there are those who scoff at the idea of making 
any present value computations on exploratory ventures. 
There are already so many uncertainties, is their 
belief, that to refine the calculation any further is 
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like tilting at windmills. The answer to this has 
been stated several times: a present value calculation 

is an estimate, but an educated one that explicitly 
recognizes time as a realistic decision consideration. 

To ignore the effect of time or to try to cope with it 
in an implict way is to incur the danger that it may 
be overlooked entirely or its effect miscalculated. 

On the basis of the concepts discussed in this section, 
the highest, lowest and most likely values of profit per barrel 
can be computed. Table 7-1 illustrates this process and is 
based on the hypothetical expenses and prices for a drilling 
venture in the Southern Louisiana on-shore area. 

The values developed in Table 7-1 will be used as the 
base figures for the sensitivity studies contained in subsequent 
Chapters . 

7 . 3 The Triangular and Beta Distributions 

It was seen in Section 7.2 that while the operator 
would probably be unable to provide the exact profit per barrel 
of oil produced he could provide his estimates of the highest 
profit per barrel he would expect, the lowest profit per barrel 
and the profit per barrel he feels most likely to occur. It 
is assumed of course that these figures will be discounted to 
reflect the time value of money. Given this data, the question 
then becomes, what distribution function best suits this set of 
parameters? If all that is known is the upper limit, the lower 
limit and most likely value (mode) it is possible to describe 
this knowledge by means of either the triangular distribution 
function or the beta distribution function. Of the two, the 



triangular function is the simplest and will be discussed first. 
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TABLE 7-1 



DISCOUNTED PROFIT PER BARREL - SOUTHERN LOUISIANA ON-SHORE 



ITEM 


MOST 

LIKELY 


HIGH 


LOW 


POSTED PRICE 


$3.25 


3.50 


3.00 


DEDUCTIONS: 








ROYALTY: 

@1/6 

@1/8 


. 54 


.44 


. 50 


SEVERANCE TAXES (see Note 1) 


. 24 


. 26 


.22 


OPERATING COSTS (including lift 
cost, labor, lease maint., 
etc . ) 


.10 


.11 


.09 


HOME OFFICE COSTS OR OVERHEAD 
(@ 5%) 


.16 


.18 


.15 


GEOPHYSICAL COSTS 


.01 


.01 


.01 


DRILLING § COMPLETION 
(see Note 2) 


.40 


.40 


.40 


LEASE ACQUISITION (see Note 3) 


. 02 


. 02 


.02 


TOTAL DEDUCTIONS 


-$1.82 


-1.62 


2.19 


NET PROFIT 


1.43 


1.88 


.81 


DISCOUNTED PROFIT (based on 
discounting at 10% for 10 
years: discount factor = 

.6144) 


$ .88 


$1.15 


$.50 


NOTE 1: [49] a. State and local severance 


taxes paid 


in 1968 



= $233,017,000. 



b. Total Petroleum liquid production in 1968 = 
2,526,008 bbl/day (average). 

c. Average field price of oil (1969) = $3.33 

$233,071,000 
2,526,008 bbl/day x 365 



= $ . 252/bbl 
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TABLE 7-1 Continued 



however, at $3.30/bbl 

Severance taxes = • 252 = 7 . 35 % 

3.30 

say 7.4% of posted price 



NOTE 2: [50] a. Assume production at 150 bbl/day. 

b. Estimated Cost of drilling and equipping wells 
in 1970 = $196,500. Allowing for inflation 

§ deeper wells assume cost of drilling = 

$217, 500/well. 

c. Assume recovery of 217,500 in 10 years then 
recovery/year = $21,750. 



$21,750 

150 bbl/day x 365 



$ . 40/bbl 



NOTE 3: a. 

b. 



c . 



d. 



Assume Lease Costs at $60/acre where a completion 
is made. 

Assume well spacing at 40 acres; 

then Lease Costs $2400/well where well is completed. 

Assume lease on unproductive land to be equivalent 
to $2400/well; 

then total lease costs = $4, 800/well. 

Assume recovery in 5 years with production rates 
= 150/bbl/well . 



— 5 — = $960/well 

$960 

150 bbl/day x 365 



.017 



say $. 02 /bbl 
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1 • 4 The Triangular Distribution Function 

The triangular distribution is defined by three values 
of the variable; i.e. the lowest valur of x = a, the highest 
value of x = c and the most likely value of x = b. The prob- 
ability distribution function w(x) can take three forms as 
shown below: 

w (b) 




Figure 7-1A Triangular Distribution - Positive Skew 




Figure 7-1B Triangular Distribution- Symmetrical 
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w(b) 




c 

Since / w(x) dx = 1 (= the area of the triangle) for a 
a 

continuous probability distribution: 

c 

/ w(x) dx = l/2w(b) (c-a) = 1 (3) 

a 

from which 



w (b) 



2 

c-a 



(4) 



In unpubilshed lecture notes Folkert Brons [11] 
proves the following properties for the triangular 
distribution : 



mean, x 



a + b + c 
3 



mode, w(b) 



c 



2 



a 



median, x 
* m 



c- 




c-a) (c-b) 
2 



(5) 

( 6 ) 

(7) 



for positive skew (see Figure 7-1A) 
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x m - a * ,/( 9 -a)(b-a) 

m V 2 



( 8 ) 



for negative skew (see Figure 7-lc) 

Brons has also derived the cumulative distribution function 
[11] and has demonstrated that: 

For x ^ b 

2 

(9) 



( x i" a ) 

W ( x < x i) = 'Cb'-a)'(c-'a) 



or , x 



^ = a + \/w(x < x^) (b-a) (c-a) 



( 10 ) 



For x < b 



(c -x • ) 

W(x > x .) = t — t-vt-I — r 
i (c-b) (c-aj 



or, x^ = c - \j~W(x > x^) (c-b) (c-a) 



( 11 ) 



( 12 ) 



Further , 



W (x < b) = 



b-a 

c-a 



(13) 



The variance for the triangular distribution is given by: 

i 2 + b 2 + c 2 - ab - ac - be 



a 2 = £ 



18 



(14) 



From which the Coefficient of Variation may be derived as 
follows : 



CV = — = 



V 



2,2 2 , , 

a + b + c - ab -ac - be 



18 

a + b + c 



x 



3 
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I 2 2 ? 

a + b + c - ab - ac - be 



CV = 



a 



+ b + c 



pv _ 'I a 2 + b 2 + c 2 - ab - ac - be 

LV rrm'Ca v 'tr+ c') 



(15) 



7 . 5 Simulation of the Triangular Distribution 

Using the data for the highest, lowest and most likely 
values of discounted profit per barrel previously developed it 
is possible to simulate the cumulative triangular distribution 
of profit per barrel W(x) as follows: 

Let a = $ . 50 
b = $.88 
c = $1.15 

x = specific values of expected discounted 
profit per barrel 
then for, $.50 .< x^ < $.88 

x^ = a + \/w(x < x^) (b-a) (c-a) 

x. = .50 + ^/w(x < x i )(. 88-.50)(1.15 -.50) 

x i = .50 + \/w (x < x~) (.2470) (16) 

for, $1.15 2, x^ $.88 

x^ = c - \/w(x > x^) (c-b) (c-a) 

Xj = 1.15 - \/w(x > x L ) (1.15-.88) (1.15-.50) 



X . 
1 



1.15 - Vw(x > x i ) (.1755) 



(17) 
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and , 

W(x < b) - |^f 

.88-. 50 

W(x £ $ . 88) = i . is- . 50 

W (x < $.88) = .585 (18) 

Further, from Equation 15 the Coefficient of Variation may be 
calculated as follows: 

^ , _ w? + + c ^ - ab - ac - be 

LV “ 1 .414 (a + b + c) (15) 



V(7 50) 2 +(.88) 2 +(1.15) 2 -(.50) (1 . 15) - ( . 88) (1 . 15) 
1.414 (.50+ .88 + 1.15) 



CV 



. 564 
3.580 



CV = .157 



A series of random numbers can now be assigned to represent 

values of W(x < x^) from .0000 to .999999 From 

Equation 16, if this random number was .585 or less the random 
number would be substituted into Equation 10, to obtain a 
value of x^. If the random number was .585 or greater 

W(x < x^) = 1.0 - W(x > x^) (19) 

would be used to obtain W(x > x^) and then that value would 
be substituted into Equation 12, to obtain the corresponding 
value of x^. This process is illustrated in Table 7-2. 
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TABLE 7-2 

ILLUSTRATION: PROFIT PER BARREL ASSUMING TRIANGULAR 

DISTRIBUTION 



Random 

No. 


W(x < x^) 


W(x > x^) 


x . = .50 + 

l 


x i = 1.15 - 


\/w(x < x^) ( . 2470) 


Vw(x > x i )(. 1755) 


.00 


.00 


(1.00) 


.50 




.01 


. 01 


(.99) 


.55 




.02 

• 


.02 

• 


(.98) 

• 


. 57 




• 

. 58 


• 

. 58 


• 

(.42) 


.87 




. 59 


(.59) 


.41 




.89 


. 60 
• 


(.60) 

• 


.40 

• 




. 89 


• 

.99 


• 

(.99) 


• 

. 01 




1— 1 

• t— 1 

i— 1 



In this way the best estimate of discounted profit per barrel, 
corresponding to a random number generated to represent a value 
of W (x < x^) , may be simulated. For use in the Monte Carlo 
Simulation, a Triangular Distribution Variable Generator was 
developed and is included in Appendix B. 

In order to test the reliability of the deviates pro- 
vided by the Triangular Distribution Generator, 1000 values of 
the profit per barrel, e, were calculated using the program 
included in Appendix B and the values of a, b, and c previously 
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discussed. The results of this test are listed on Table 7-3 
under the heading of Triangular Distribution. The values of 
profit per barrel thus simulated were plotted on Figure 7-2 as 
circles. Also included on Figure 7-2, as triangular points, 
are values of profit per barrel calculated manually using the 
technique outlined on Table 7-2. It is apparent from Figure 
7-2 that there is excellent correlation between the manually 
calculated values and those calculated using the Triangular 
Distribution Generator, thus confirming the reliability of the 
Generator as a simulation tool. Dispite this excellent correl- 
ation one item should be recognized. By the very nature of the 
random number generator which generates values of numbers used 
to simulate W (x < x^) a value of 1001 for W (x < x^) can never 
be achieved. Instead the highest value of W(x < x^) we can 

achieve using the generator is .9999 Thus, there is no 

possibility that the highest value of profit per barrel (e) 
can ever be reached. This shortcoming is considered minor 
however, considering the excellent correlation for virtually 
the entire range of deviates. 

7 . 6 The Beta Distribution Function 

In this derivation of the Beta distribution function 
Brons [13] notes that a disadvantage of the Triangular distri- 
bution is that its derivate ^ ^ x — becomes undetermined at x = 

dx 

a, b and c, i.e. the only points where some information about 
x is available. This means that at the outer limits a and b, 




.50 .60 .70 .80 .90 1.00 1.10 

Discounted Profit/Barrel ($) 



100 



Figure 7-2 

Cumulative Distribution Curve - Triangular Distribution 



rirtbtitb 
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w(x) suddenly drops to zero. This is depicted in Figure 7-3. 




i 





1 

dw (x) , 

dx 

i 







x 



Figure 7-3 - Triangular Distribution Derivative 

He also notes that a more gentle dropping off of the probabil- 
ity of occurence of x near the end points often better ex- 
presses reality - in that occurances in a real situation would 
normally taper off rather than end abruptly. This concept can 
be expressed by the Beta Distribution shown in Figure 7-4. 
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w(x) 



a 



b 



c 



■»— x 



Figure 7-4 Beta Distribution Function 



The Beta distribution, like the Triangular distribution, is 
defined by the three parameters a, b and c. 

The Beta Distribution is given by 



of gamma functions where a and g are defined by the values of 
a, b, and c. A derivation of the relationships between a and 
g and a, b and c (and of the Beta distribution itself) is be- 
yond the scope of this thesis and can be found in several 
excellent discussions of this distribution [13] [19] . 

In an unpublished report [30], W.G. Lesso developed a 
computer program which generates deviates having a Beta dis- 
tribution. The program uses as inputs, the mean and variance 
of x (as expressed in a PERT type analysis) : 




( 20 ) 



It can be seen that the Beta distribution is a ratio 
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.3. _ a + 4b + c 


6 


(21) 


a 2 = I _ c _-. a .). 2 


3C 36 


(22) 


From which the mean and variance 
bution may be obtained as follows: 


for the Beta distri- 


! 1 

X • 
1 

p 


x b c-a 

2 °x 2 


(23) 


b f ,2 

(c-a) 


(24) 



The Lesso developed program, modified for use in this study, 
is included in Appendix C. 

For comparision purposes the Beta Distribution Variable 
Generator was used to calculate 1000 values of the profit per 
barrel, e, using the same parameters of a,b and c as were used 
to calculate the profit per barrel using the Triangular Dis- 
tribution Variable Generator. The values of profit per barrel 
thus calculated are listed in Table 7-3 under the heading of 
Beta Distribution. Cumulative frequency curves for the Beta 
and Triangular distributions, using the data contained in Table 
7-3, is plotted on Figure 7-5. It can be seen from Figure 7-5 
that very little difference exists between the distributions 
except for values at the lower end of the curves. It should be 
noted, however, that the mean and standard deviation for the 
Beta distribution was computed using the values as expressed in 
a PERT type analysis (see Equations 21 and 22) rather than in 
the form used in the Triangular distribution (see Equations 5 
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and 14) . 



Because of its relative simplicity a Triangular 
distribution will be assumed for the variable, profit per 
barrel. In the sensitivity studies which follow however, the 
Beta distribution will be substituted for the Triangular dis- 
tribution and its effect on the venture profitability will be 
analyzed. 
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TABLE 7-3 

COMPUTER SIMULATED DEVIATES FOR DISCOUNTED PROFIT PER BARREL 
ASSUMING TRIANGULAR AND BETA DISTRIBUTIONS 



parameters 

lowest profit per barrel: $.50 
most likely profit per barrel: $.88 
highest profit per barrel: $1.15 



rof it/bbl 
($) 


TRIANGULAR 


DISTRIBUTION 


BETA DISTRIBUTION 


Relative 

Frequency 

w = £ - - 


Cumulative 

Frequency 

W(x < x^) 

(%) 


Relative 

Frequency 

f 


Cumulative 

Frequency 

W (x < Xj,) 

m 


1000 

(%) 


Tooo 

m 


50 


0 


0 


0 


0 


51 


0 


0.0 


0 


0.0 


52 


0.2 


0.2 


0 


0.0 


53 


0 


0.2 


0 


0.0 


54 


0.5 


0.7 


0 


0.0 


55 


0.9 


1.6 


0 


0.0 


56 


0.6 


2.2 


0.1 


0.1 


57 


0.6 


2.8 


0 


0.1 


58 


0.7 


3.5 


0 


0.1 


59 


0.7 


4.2 


0.1 


0.2 


60 


0.7 


4.9 


0.3 


0. 5 


61 


0.7 


5.6 


0 


0.5 


62 


0.6 


6.2 


0.6 


1.1 


63 


1.1 


7.3 


0.6 


1.7 


64 


1.3 


8.6 


0.6 


2.3 


65 


1.4 


10.0 


0.8 


3.1 


66 


1.0 


11.0 


1.0 


4.1 


67 


1.3 


12.3 


0.9 


5.0 


68 


1.2 


13.5 


1.0 


6.0 


69 


0.8 


14.3 


1.4 


7.4 


70 


1.9 


16.2 


1.2 


8.6 


71 


1.9 


18.1 


2.0 


10.6 


72 


1 . 4 


19.5 


2.2 


12.8 


73 


1.7 


21.2 


2.0 


14.8 


74 


1 . 9 


23.1 


1.6 


16.4 


75 


1.9 


25.0 


1.8 


18. 2 


76 


2.0 


27.0 


2.6 


20.8 
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TABLE 7-3 Continued 



Prof it/bbl 
($) 


TRIANGULAR 


DISTRIBUTION 


BETA DISTRIBUTION 


Relative 

Frequency 

w - f 


Cumulative 

Frequency 

W(x < x i ) 

m 


Relative 

Frequency 


Cumulative 

Frequency 

W (x < x^ 

(%) 


1000 

(*) 


w TTJ0TT 

m 


77 


1.8 


28.8 


2.9 


23.7 


78 


2.5 


31.3 


1.6 


25.3 


79 


2.5 


33.8 


3.4 


28.7 


80 


2.1 


35.9 


2.3 


31. 0 


81 


2.3 


38.3 


2.8 


33.8 


82 


2.5 


40.7 


3.1 


36.9 


83 


3.0 


43.7 


2.3 


39.2 


84 


3.2 


46.9 


3.2 


42.4 


85 


2.6 


49.3 


3.1 


45.5 


86 


2.8 


52.3 


3.4 


48.9 


87 


3.2 


55.5 


3.7 


52.6 


88 


3.9 


59.4 


3.6 


56.2 


89 


2.6 


62.0 


3.3 


59.5 


90 


2.9 


64.9 


4.3 


63.8 


91 


2.6 


67.5 


3.7 


67.5 


92 


1.9 


69.4 


3.2 


70.7 


93 


3.2 


72.6 


2.1 


72.8 


94 


2.5 


75.1 


3.5 


76.3 


95 


2.5 


77.6 


2.1 


78.4 


96 


1.8 


79.4 


3.3 


81.7 


97 


2.5 


81.9 


2.3 


84.0 


98 


1.7 


83.6 


3.1 


87.1 


99 


2.2 


85.8 


1.6 


88.7 


100 


2.2 


88.0 


2.1 


90.8 


101 


1.5 


89.5 


1.4 


92.2 


102 


2.0 


91.5 


1.5 


93.7 


103 


1.3 


92.8 


1.3 


95.0 


104 


1.4 


94.2 


1.2 


96.2 


105 


0.8 


95.0 


0.6 


96.8 


106 


1.5 


96.5 


0.7 


97.5 


107 


0.6 


97.1 


1.0 


98 . 5 


108 


0.9 


98.0 


0.5 


99.0 


109 


0.4 


98 . 4 


0.3 


99.3 


110 


0.5 


98.9 


0.2 


99.5 
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TABLE 7-3 Continued 



Prof it/bbl 
($) 


TRIANGULAR 


DISTRIBUTION 


BETA DISTRIBUTION 


Relative 

Frequency 

w = f 


Cumulative 

Frequency 

W(x < x i ) 

(*) 


Relative 

Frequency 

w = £ 


Cumulative 

Frequency 

W(x <_ x^) 

(%) 


1000 

m 


1000 

(%) 


111 


0.2 


99.1 


0.3 


99.8 


112 


0.5 


99.6 


0 


99.8 


113 


0.4 


100.0 


0.2 


100.0 


114 


0 


100.0 


0 


100.0 


115 


0 


100.0 


0 


100.0 



Figure 7-5 

Triangular and Beta Distribution Curves 
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CHAPTER VIII 



DRY HOLE COSTS 

In Chapter V we discussed the concept of probability 
of success and how it related to the number of successful wells 
to be expected in a drilling program. To account for the num- 
ber of wells that are unsuccessful we have included dry hole, 
or unsuccessful well, costs into our profitability model. 

While it can be argued that dry hole costs should have 
a continuous frequency distribution, just as the other parameters 
in the model, we will assume for simplicity sake that the dry 
hole costs are constant in this thesis. 

To be consistent with the values developed in Table 
7-1, we will assume that the drilling venture is to take place 
in the Southern Louisiana on-shore area. From January to Dec- 
ember 1970 there were a total of 1245 wells drilled in this 
area with a total footage of 12,703,309 [51]. Therefore, the 
average well depth would be equal to 10,200 feet/well. Using 
1969 average dry hole costs of $13. 23/ft [50], the average dry 
hole cost per well would be $135,000. Allowing for inflation 
and deeper wells we will assume a constant dry hole cost of 
$150,000 for each unsuccessful well in the program. 
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CHAPTER IX 



COMPUTER PROGRAM 

The FORTRAN IV computer program, as contained herein, 
was developed to apply the Monte Carlo simulation technique to 
an exploratory well drilling venture using the concepts and 
relationships established in previous chapters, to determine 
the best estimate of discounted profit and to accumulate data 
for each iteration so that a cumulative probability curve can 
be plotted. 

The program' is designed to call up the same series of 
random numbers with each separate computer run so that sen- 
sitivity studies can be made. The random number generator 
used in this program is a function designed for the CDC 6600 
computer system and is called by RANF(O). The generator pro- 
duces pseudo random numbers, distributed uniformly in the in- 
terval 0<x<l. 

The program accepts as data: the number of wells 

drilled, probability of success, dry hole costs, and the value 
M, which is the number of end-points entered in Line 5 plus 
one, (Line 4); the median value of Q, the value of Q at the 
median + 1 standard deviation, Q84, and the boundary conditions, 
Q max. and Qmin., (Line 6); and the lowest, most likely and 
highest profit per barrel, A, B and C, (Line 7). In addition, 
the end points of the random number ranges corresponding to the 
probability of the number of successful wells is entered as NEP , 
(Line 5) . 
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With this data entered, values are initialized (Lines 
12, 13, 14, 15, 18, 20) and the Standard Deviation and Mean of 
the Lognormal distribution are calculated (Lines 16 and 17) . 

For each iteration specified (Line 23) the first series 
of operations establishes the number of successful wells (Lines 
24-32). To do this a random number is called, RANF(O), mult- 
iplied by 10,000 and compared with the end points corresponding 

to the number of successes previously entered (Line 5). Thus, 

* 

the random number called establishes the number of successes, 
NSW. 

If the number of successes is greater than 0, the pro- 
gram is designed to calculate the size reserve for each success 
and to total them to obtain the value QSUM, which is the total 
reserve size discovered (Lines 37-49). For a detailed descrip- 
tion of the Lognormal distribution generator see Appendix A. 

With the total reserve size determined the profit per 
barrel, for this iteration, E, assuming a triangular distribu- 
tion, is calculated (Lines 51-57). For a detailed discription 
of the Triangular distribution ganerator see Appendix B. 

Having calculated or inputed all the parameters the 
iteration profit, PROFIT, is computed (Line 58) using the re- 
lationship expressed in Chapter III, Equation 4. Lines 59 and 

60 accumulate this iteration profit and average the value, PTOT, 
over the number of iterations run, MM. XSQ and XSQTOT in Lines 

61 and 62 are used to calculate the standard deviation and coe- 
fficient of variation later in the program. 
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Lines 64 through 176 are inserted to store the values 
of iteration profit in arrays by PROFIT intervals in order to 
compute the probability of the interval, PROB, and the cumu- 
lative probability, CUMP (Lines 189-202). 

The variance, VAR, standard deviation, SD, and coeff- 
icient of variation, CV, are calculated for the particular set 
of model parameters in Lines 178, 179 and 180 using information 
previously determined in Lines 61 and 62. 

A listing of the FORTRAN IV computer program is con- 
tained on the following pages. 





PROGRAM PROFIT ( INPUT .OUTPUT ) 


1 




DIMENSION N£P (25) ♦ IY(40). QLOGN(20>» D(40). YY(40>. PROB<40>. 


2 




1 CUMP (40 ) 


3 




READ 100. N.P.CU.M 


4 




READ 110. (NEP(K), K=1,M) 


5 




READ 111. QMEDIAN. Q84, UMAX. QMIN 


6 




READ 112. A, B. C 


7 




PRINT 155 


0 




PRINT 113 


9 




PRINT 1 1 A . N. P. CU. QMEDIAN. Q84. QMaX. QMIN 


10 




PRINT 160. A. B. C 


11 




PROF = 0.0 


12 


20 


PROFIT = 0.0 


13 




XSQ = 0.0 


14 




XSQTOT = 0.0 


15 




SO = ALOG (Q84) - ALOG (QMEDIAN) 


16 




TMU = ALOG (QMEDIAN) 


17 




MM = 1 


10 




PRINT 115 


19 




DO 300 I = 1.40 


20 




I Y ( I ) = 0.0 


21 


300 


CONTINUE 


22 




DO 99 II = 1.1000 


23 




CALL RANF(O) 


24 




Wl = RANF(O) 


25 




W = Wl«10000. 


26 




J = W 


27 




DO 10 I = 1 * M 


28 




IF < J.Le.N£P( I) ) GO TO 11 


29 


10 


CONTINUE 


30 


11 


NSW = 1-1 


31 




IF (NSW. FQ. 0.0) 32.40 


32 


32 


QSUM = 0.0 


33 




E = 0.0 


34 




GO TO 56 


35 


40 


CONTINUE 


36 




QSUM = 0.0 


37 




DO 45 I = 1 .NSW 


30 




RA = RANF(O) 


39 




RB = RANF(O) 


40 




V = (-2.0* ALOG (RA))**. 5* C0S(2. 0*3.141 6* RB) 


41 




QNORM = V*SD ♦ TMU 


42 




QLOGN(I) = EXP (QNORM) 


43 




IF (QLOGN(I)-QMIN) 46.47.48 


44 


46 


QLOGN(I) = QMIN 


45 




GO TO 47 


46 


40 


IF (QLOGN(I)-QMAX) 47.47.49 


47 


49 


QLOGN(I) s QMAX 


48 


47 


QSUM = QSUM ♦ QLOGN(I) 


49 


45 


CONTINUE 


50 




CALL RANF(O) 


51 




WXLTXI = RANF(O) 


52 




WXGTXI = 1.0 - WXLTXI 


53 




XI = A ♦ SQRT(WXLTXI*(B-A)*(C-A) ) 


54 




IF (XI.LE.B) 55.50 


55 


50 


XI = C - SQRT(WXGTXI*(C-B)*(C-A) ) 


56 


55 


E = XI/100. 


57 


56 


PROFIT = QSUM*E - <N-NSW)*CU 







PhOF = PROF I r ♦ PROF 






59 




PTOT = P&OF/MM 








60 




XSQ = ( PROF I T ) *■*2 








61 




XSQTOT = XSQTOT ♦ XSQ 






62 




PRINT life, MM, NSW, 


QSUM, 


Ft PROFIT, PTOT 


63 




MM = MM ♦ I 








64 




IF (PROFIT ,LT ,-10,0) 


GO 


TO 


61 


6b 




I F ( PROF I T . L T . 0.00) 


GO 


TO 


62 


66 




IF (PROFIT.LT. 10.0) 


GO 


TO 


6 3 


67 




IF (PROFIT. LT. 20.0) 


GO 


TO 


64 


68 




IF (PROFIT.LT. 30.0) 


GO 


TO 


65 


69 




IF (PROFIT. LT. 40.0) 


GO 


TO 


66 


70 




IF (PROFIT.LT. 50.0) 


GO 


TO 


67 


71 




IF (PROFIT.LT. 60.0) 


GO 


TO 


68 


7? 




IF (PROF IT .LT. 70.0) 


GO 


TO 


69 


73 




IF (PROFIT.LT. 80.0) 


60 


TO 


70 


74 




IF (PROFIT.LT. 90.0) 


GO 


TO 


71 


7b 




IF (PROFIT.LT. 100.0) 


GO 


TO 


7? 


76 




if (PROfit.lt. no. 0) 


GO 


TO 


73 


77 




IF (PKOFIT.LT. I 20.0) 


GO 


TO 


74 


78 




IF (PROF IT.LT. 130.0) 


f,0 


TO 


75 


79 




IF (PROF IT .LT. 140.0) 


GO 


TO 


76 


80 




IF (PROF IT .LT . 150 .0) 


GO 


T 0 


77 


81 




IF (PROFIT.LT. 160.0) 


GO 


TO 


78 


82 




IF (PROFIT.LT. 170.0) 


GO 


TO 


79 


83 




IF(PROFIT.LT. 180.0) 


GO 


TO 


80 


84 




IF (PROFIT .LT. 190.0) 


GO 


TO 


8] 


85 




IF (PROFIT .LT. 200.0) 


GO 


TO 


82 


86 




IF (PROFIT .LT. 210.0) 


GO 


10 


83 


8/ 




IF (PROF IT .LT. 220.0) 


GO 


TO 


84 


88 




IF (PROF IT .LT. 230.0) 


r;0 


TO 


85 


69 




IF (PROF IT .LT .240.0) 


GO 


TO 


86 


90 




IF(PROFIT.LT. 250.0) 


GO 


TO 


87 


91 




IF (PROFIT.LT. 260.0) 


GO 


TO 


88 


92 




IF (PROF IT.LT. 270.0) 


GO 


TO 


89 


93 




IF (PROFIT.LT. 280.0) 


GO 


TO 


90 


94 




IF (PROF IT.LT. 290.0) 


GO 


TO 


91 


95 




IF (PROFIT .LT. 300.0) 


GO 


TO 


9- 


96 




IF (PROFIT.LT. 350.0) 


GO 


TO 


93 


97 




IF(PR0FIT.LT. 400.0) 


GO 


TO 


94 


98 




IF (PROF IT .LT. 450.0) 


GO 


TO 


95 


99 




IF (PROFIT .LT .500.0) 


GO 


TO 


96 


100 




IF (PROFIT.LT. 550.0) 


GO 


TO 


97 


101 




GO TO 99 








102 


61 


I Y ( I ) = I Y ( 1 ) ♦ I 








103 




GO TO 99 








104 


62 


I Y ( 2 ) = I Y ( 2 ) ♦ I 








105 




GO TO 99 








106 


63 


I Y (3) = I Y (3) ♦ I 








107 




GO TO 99 








I 08 


64 


I Y (4) = IY (4) ♦ I 








109 




GC TO 99 








I 10 


65 


I Y ( 5 ) = I Y ( 5 ) ♦ I 








111 




GO TO 99 








I 12 


66 


IY (6) = I Y (6) ♦ I 








113 




GO TO 99 








114 


67 


I Y ( 7 ) = I Y ( 7 ) ♦ 1 








115 




GC TO 99 








116 



68 

69 

70 

71 

72 

73 

7m 

7b 

76 

7 7 

78 

79 

80 

81 

82 

83 

88 

85 

66 

87 

88 

89 

90 

91 

92 

93 

98 

95 

96 



1 1 7 

118 

1 19 

120 

121 

122 

123 

128 

125 

126 

127 

128 

129 

130 

131 

132 

133 

1 38 

135 

136 

137 

138 

139 

180 

181 

182 

183 

188 

185 

186 

187 

188 

189 

150 

151 

152 

153 

158 

155 

156 

157 

158 

159 

160 

161 

162 

163 

168 

165 

166 

167 

168 

169 

170 

171 

172 

173 

178 



I Y (8) 


= 


I Y ( 8 ) ♦ 


1 




uO ro 


99 








I r (9) 


- 


I Y ( 9 ) ♦ 


1 




GC TO 


99 








IY ( 10 ) 


= 


IY ( 10) 


♦ 


1 


go ro 


99 








1Y ( 11 ) 


- 


IY(ll) 


♦ 


1 


GO 1 0 


99 








1Y ( 12) 


= 


1Y( 12) 


♦ 


1 


GO 10 


99 








1 Y ( 13) 


- 


I Y ( 13) 


♦ 


1 


GO TO 


99 








I Y ( 18) 


= 


I Y ( 18) 


♦ 


1 


GC TO 


99 








1 Y ( 15) 


- 


I Y ( 15) 


♦ 


1 


TO 


98 








I Y ( 16) 


= 


I Y ( 16) 


♦ 


1 


r o 


99 








1 Y (1 7) 


= 


I Y ( 17) 


♦ 


1 


uO 1 0 


99 








1 Y ( 18) 


= 


I Y ( 18) 


♦ 


1 


mO 10 


99 








I Y ( 19) 


= 


I Y ( 19) 


♦ 


1 


GC TO 


99 








IY (20) 


= 


1 Y (20) 


♦ 


1 


GC TO 


99 








I Y (21 ) 


= 


1 Y (21) 


♦ 


1 


,C 10 


^9 








1 Y (22) 


= 


IY (22) 


♦ 


1 


<jO TO 


99 








It (23) 


=. 


IY (23) 


♦ 


1 


(jrC 10 


99 








I Y (28) 


- 


I Y (28) 


♦ 


1 


GC TO 


99 








IY (25) 


= 


I Y (25) 


♦ 


1 


GC TO 


99 








I Y (26) 


= 


IY (26) 


♦ 


1 


bO TO 


99 








I Y (27) 


= 


I Y (27) 


♦ 


1 


Go TO 


99 








1 Y (28) 


= 


IY (28) 


♦ 


1 


GC TO 


99 








1 Y (29) 


= 


1 Y (29) 


♦ 


1 


GO TO 


99 








I Y (30) 


= 


I Y (30) 


♦ 


1 


GO TO 


99 








IY ( 31 ) 


= 


IY (31 ) 


♦ 


1 


GO TO 


99 








I Y ( 32) 


= 


IY (32) 


♦ 


1 


GO TO 


99 








IY (33) 


= 


IY (33) 


♦ 


1 


GC TO 


99 








I Y (38) 


= 


IY (38) 


♦ 


1 


GC TO 


99 








IY (35) 


= 


1 Y ( 35 ) 


♦ 


1 


GC TO 


99 








IY (36) 


= 


IY (36) 


♦ 


1 


GO TO 


99 
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99 

60 



210 



£20 

230 

200 

100 

no 

in 

112 

113 



114 

115 



116 

140 



141 

144 

145 

150 

155 

160 



1Y<37) = I Y ( 37 ) ♦ 1 
CONTINUE 

MM = MM - 1 

VAR = ( XSQTOT - (MM*(PT0T**2) ) )/(MM - 1) 
SO = SORT (VAR) 

CV = SD/PTOT 



PR 


1NT 


140, 


PR 


INT 


141 , 


PR 


I NT 


144 


PR 


INT 


155 


PR 


INT 


113 


PR 


INT 


114, 


PR 


INT 


160, 


PR 


INT 


145 


SUMM = 


0.0 


D ( 


1) = 


-10 


DO 


200 


IM 



MM * PTOT 
VAR* SD * CV 



N , P, CU, QMEOIAN* Q84, UMAX , QM I N 
A * B , C 



0 

1*37 



YY ( IM) = I Y ( IM) 

PROS ( IM) = YY(IM)/1000. 
SUMM = SUMM ♦ PRO^(IM) 
CUMP(IM) = SUMM 



Kj = IM ♦ 1 

lF(IM-32) 210, 220 , 220 
D (KJ) = 0( IM) ♦ 10.0 
GO TO 230 

D (KJ) = D( IM) ♦ 50.0 
CONTINUE 

PRINT 150* D(IM), PROp(IM)* CUMP(IM) 

CONTINUE 

FORMAT ( 10X, 15, 2F 10.3, 15) 

FORMAT (10Xfl7l4) 

FORMAT ( 1 OX , 4F 1 0 • 3 ) 

FORMAT (10X,3F10.2> 

FORMAT ( /// , 6X , *N0 • OF W£LLS*,3X, *PR0B 0F*,2X,*0RY HOLE* , / , 9X , *DR I 
1LLED*,5X,*SUCCESS*,4X,*C0ST*,4X,*GMEDIAN*,3X,*G84*,5X,*QMAX*,5X, 
2*QM IN* ) 

FORMAT ( 1 OX , I5,8X,F4.3,4X,F4.3 , 3X , F7 . 2 , 2X , F 7 . 2 , 2X , F 7 . 2 , 2X , F7 . 2// ) 
FORMAT <5X,* ITERATION*, 3X,*NUMBER OF* , 4X , *F IELO* , 2X , *PROF I T PER*,2X 
1 ** ITERATION* ,3X,*BEST ES T I M A TE* , / 1 7X , *SUCCESSES* , 4X , *S I ZES* * 5X t *B A 
2RREL*,6X,*PR0 f1T*,6X,*0f PROF IT*,//) 

FORMAT ( 5X, 14,1 IX, 13 ,6X,F6.2,6X,F4.2*7X,F6.2f8X,F6.2) 

FORMAT <///,5X,*THE BEST ESTIMATE OF PROFIT BY THE MONTE CARLO SIM 
1ULATION METHOD**/7X,*AT *,I4(* ITERATIONS IS *,F7.2,* MILLION DOLL 
1 AkS* ) 

FORMAT ( /// ,5X , *THE VARIANCE IS * , F 1 2 . 3 , / ,5X , *THE STANDARO OEVIATI 
1 ON IS *,F12.3*/*5X,*THE COEFFICIENT OF VARIATION IS *,F6.4> 

FORMAT ( 1 H 1 ) 

FORMAT ( 1 5X , *UPPER LIMIT OF* , 7X , *PROB AB IL I T Y* * 5X , *CUMMULAT I VE* , 
1 / , 1 5X , *PROF I T INTERVAL*, 5X,*0F INTERVAL* ,5X,*PR0BABILITY*,///) 
FORMAT (^0X,F4.0 t 15X,F5.4tl 1X,F5.4) 

FORMAT (1H1,15X,*PR0FIT DISTRIBUTION OF AN EXPLORATORY DRILLING PR 
10GRAM*) 

FORMAT ( / 1 5X , * THE TRIANGULAR DISTRIBUTION FOR DETERMINATION OF THE 
1*,/15X,*PR0FIT PER BARREL HAS THE FOLLOWING EXTREMES IN CENTS - *, 
2/ * 1 5X , * A = *,F6.2,7X,*B = *,F6.2,7X,*C = *,F6.2,///) 

End 
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CHAPTER X 



BEST ESTIMATE OF DISCOUNTED EXPECTED PROFIT USING BASE DATA 



In previous chapters the profitability model was dev- 
eloped and the individual parameters and their frequency dis- 
tributions were discussed. In addition, certain Base Data was 
established for each parameter so that the best estimate of 
discounted expected profit, E(P), could be calculated. This 
Base Data is summarized in Table 10-1. 



TABLE 10-1 
BASE DATA 



Parameter 


Distribution 


Value 


Number of Wells Drilled - n 


- 


20 


Probability of Success - p 


Binomial 


.15 


Reservoir Size 


Lognormal 




Q at 51 probability 




5 x 10 6 bbl 


Q at 951 probability 




79 x 10 6 bbl 


Q maximum 




140 x 10 6 bbl 


Q minimum 




1 x 10 6 bbl 


Discounted Profit/Barrel 


Triangular 




Lowest - e„ 

a 




$ . 50/bbl 


Most Likely - e^ 




$ . 88/bbl 


Highest - e c 




$ 1 . 15/bbl 


Dry Hole Cost - C u 


None 


$ 150,000 


Number of Iterations 




1,000 



84 



85 



On the basis of this Base Data and using the Fortran 
IV program included in Chapter IX a value of Best Estimate of 
Discounted Expected Profit, E(P), of $68.02 million was obtained. 
See Appendix D for a printout of the number of successes, 
reserve sizes, profit per barrel, iteration profit and average 
best estimate of discounted profit for each of the 1000 iter- 
ations using the Base Data. In addition, Appendix D contains 
cumulative probability data for each specified interval of E(P). 
Using this cumulative probability information the histogram 
included as Figure 10-1 was prepared from which it can be seen 
that : 

1. The Median Value (at 50% probability) of Discounted 
Profit is $56.03 million. 

2. With 90% confidence we could expect a discounted 
profit of between $1.51 million and $170.0 million. 

3. The probability of a loss is 4.2%, or stating it 
another way, there is a 95.8% chance of a profit. 

Further, from the printout in Appendix D we see that: 

4. The standard deviation is $53.66 million. 

5. The coefficient of variation is .7889. 

The significance of the standard deviation and coefficient of 
variation as measures of risk were discussed in Chapter IV. 



866 



Figure 10-1 

Cumulative Probability - Base Data 
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0 



40 



240 



280 



99.8 99.9 




CHAPTER XI 



SENSITIVITY ANALYSIS 

Freedom from uncertainty is a luxury rarely enjoyed by 
the contemporary decision maker. In a society characterized by 
change, uncertainty is an established and accepted fact of life. 
This is especially true in the oil business. Notwithstanding 
the difficulties, the decision maker is expected to attain 
satisfactory results as measured by the goals of his organiza- 
tion . 

In the face of uncertainty, the most recurring quest- 
ions he faces are of the form, "What if...? 

1. What if labor costs rise and dry holes cost twice 
as much? Three times as much? 

2. What if the posted price of oil decreases by 301? 

3. What if we drill 30 exploratory wells rather than 
20? How about 10 wells? 

The "what if" question may be viewed as an introduction to ths 
sensitivity analysis. In its mathematical sense, sensitivity 
analysis is a study to determine how possible changes or errors 
in parameter values affect model outputs. In its applied or- 
ganizational setting, sensitivity analysis may be broadly de- 
fined as a study to determine the responsiveness of the con- 
clusions of an analysis to changes or errors in parameter 
values used in the analysis. In other words, would our decision 
be any different if a parameter value, such as the probability 
of success in out profitability model, were to change. The 



87 



88 



prime function of sensitivity analysis can be said to be to 

provide a better understanding of risk. Sensitivity analysis 

can provide guidelines for allocating organizational resources 

to data collection and data refinement activities [40] . The 

benefits of subjecting models to sensitivity analysis are 

summarized by Arnoff and Netzorg [2] : 

The use of an operations research model is especially 
important and advantageous in that it permits experi- 
mentation "on paper", without manipulation of the act- 
ual system. In using the model, one can assess the 
sensitivity (response) of the system to a wide variety 
of conditions - without requiring either the time, 
expense, or risks associated with experimenting with 
the system itself (if such experimentation would, in 
fact, be possible and meaningful). Hidden relation- 
ships can be brought to light and brought to bear upon 
decisions and control of activity. 



In this thesis a profitability model for an exploratory 
well drilling venture was developed and tested with Base para- 
meters. A decision problem is now faced. If the model was 
developed in order to provide results which would serve as a 
basis for recommending action to management, it is possible, 
at this point, to use the model results as is, or to select 
other alternatives: to collect additional data; refine initial 

parameters: or even alter the model itself. On the significance 
of this decision, Willian T. Morris [35] comments: 

The nature of management science is such that sooner or 
later it has to recommend some action other than making 
more observations. Usually this comes sooner rather 
than later. Eventually it will suggest that action be 
taken consistent with either acceptence or rejection 
of the hypothesis. It will recommend management "act 
as if" the hypothesis were true or "as if" it were 
false. Indeed, the fundamental question in management 
science may well be when to stop collecting data and 
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developing the model, and when to produce a recommen- 
dation for action. 

In the past the emphasis has been on whether information re- 
quired for decision making could be developed. With our grow- 
ing capacity to generate information or data, the emphasis now 
appears to be shifting to the question as to what extent this 
capacity can be utilized beneficially. It may be recognized 
that while the collection and refinement of data are costly, 
it may serve to improve recommentations for action. In most 
cases, however, little basis exists for projecting whether the 
additional information is worth the cost. Statistical decision 
theory, and a model such as the one developed in this thesis, 
will assist in the decision as to whether additional informa- 
tion is "worth the cost". The basic requirements for the stat- 
istical decision theory model are summarized by Morris [36] as 
follows : 

[The Model] must be able to enumerate the possible 
outcomes of future data collection efforts and, further, 
to compute the probabilities of these outcomes. In 
addition, it must be possible to indicate just how the 
information willquantitatively change the decision 
maker's view of his choice. 

In an excellent example, which illustrates this concept, C. 
Jackson Grayson [21] , demonstrates the application of decision 
theory to drilling decisions by oil and gas operators and also 
to the decision problem of whether to purchase additional in- 
formation before making a final drilling decision. 

Clearly then, a sensitiviey analysis which discloses 
decision insensitivity to estimated parameter variations 
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precludes the necessity of introducing more costly, time 
consuming data, collection data refinement and statistical 
decision models into the analysis. 

It is in this vein that the sensitiviey studies are 
conducted in this thesis. The profitability model provides 
the necessary vehicle for the analysis. From the results it 
is possible to determine the outcomes, as well as probabilities 
of these outcomes, associated with changes in parameter values 
thus providing some answers to the "What if?" questions posed 
at the beginning of this chapter. 



CHAPTER XII 



SENSITIVITY STUDIES 
12.1 Outline of Studies 

The sensitivity studies were conducted using the 
profitability model developed in Chapter III. In these studies 
each parameter in the model was varied from the base values as 
summarized in Table 10-1. Each time a parameter was varied all 
remaining parameters were set to their base values and the FOR- 
TRAN IV program, developed in Chapter IX, was executed 1000 
times. In this manner a printout similar to the one for the 
Base data, included in Appendix D, was obtained for each para- 
meter change or Case. 

Thus, for each Case, at 1000 iterations, we have the 
best estimate of discounted profit, the standard deviation, the 
coefficient of variation, the iteration profit, cumulative 
probability data, etc. Tables 12-1 through 12-9, grouped by 
parameter being varied, contain a summary of some of this data 
Cat 1000 iterations) as well as extrapolated information such 
as the median discounted profit, the maximum and minimum dis- 
counted profits within 90% confidence limits and the chance 
of making a profit greater than zero. The values using the 
Base data have been included on each Table to assist in eval- 
uation. It should be noted that each time the model was run 
the same series of random numbers were called up so that com- 
parisons between Cases could be made. Figures 12-1 through 
12-8 contain histograms for several Cases, plotting the cum- 
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ulative probability for each specified interval of discounted 
expected profit. Figures 12-1 through 12-8 are similar in 
presentation to Figure 10-1 for the Base data. On each hist- 
ogram the 901 Confidence Limits have been represented by dashed 
lines at 5% and 95% cumulative probabilities. 

To assist in the comparison of each Case to the Base, 
Tables 12-1A through 12-9A have been prepared and are inserted 
behind the Table it serves to amplify. These Tables contain 
the percentage differences between the Case parameter set 
values and the Base parameter set values. The percentages 
were calculated using the following relationship: 



A (in %) 



(Case data value) - (Base data v alue - ) 

(Base data value) x 100, (1) 



In addition to the Cases noted above, the Base data 
was run for 10,000 iterations rather than 1000 iterations. 
Additionally, the Beta Distribution Variable Generator was 
substituted for the triangular distribution of profit/barrel 
in the FORTRAN IV Program and the Base data was run for 1000 

iterations. The results of both these variations are summar- 
ized on Table 12-10. 



12.2 General Observations 

One of the prime purposes of the sensitivity studies 
is to assess the change in observed values resulting from a 
change in parameter value. The majority of the effects of 
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Figure 12-1 

Cumulative Probability - ’n' Parameter Set 
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Figure 12-2 

Cumulative Probability p 1 Parameter Set 
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Figure 12-3 

Cumulative Probability - Q at 51 Parameter Set 
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Figure 12-4 

Cumulative Probability - Q at 951 Parameter Set 



115 



x 

Ad 

X 

s 




o> 

o» 



IT) 

<T> 



O 

o> 



O 

k*> 



O 

cn 



o 

M 



«n 



CM 



o 40 80 120 160 200 240 280 

PROFIT (In Million Dollars) 



x 

vH 

X 



99.8 99.9 



99.8 



116 



Figure 12-5 

Cumulative Probability - C u Parameter Set 
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Figure 12-6 

Cumulative Probability - e(a) Parameter Set 
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Figure 12-7 

Cumulative Probability - e(b) Parameter Set 
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Figure 12-8 

Cumulative Probability - e(c) Parameter Set 
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these parameter changes can be observed on Tables 12-1 through 
12-9 and 12-1A through 12-9A and require no further datailed 
discussion. The analysis of a parameter change to a change in 
Best Estimate of Discounted Profit and to the risk as measured 
by the changes in Coefficient of Variation will be discussed 
in the next Section. 

At this point however, some general observations 
regarding the sensitivity of the model to parameter changes 
can be made. First, an inspection or the Tables will reveal 
that the one single observed value which fluctuated the great- 
est on a percentage basis was the minimum value of discounted 
profit which could be expected with 90% confidence. The reason 
for this high percentage fluctuation was the relatively small 
numbers involved. Thus, a variation of $7.87 million on a 
$1.51 million Base resulted in a percentage change of approx- 
imately 980%. 

The second observetion made from the Tables is that, 
except in isolated cases, the chance of a profit greater than 
zero was fairly constant at approximately 95 to 96%. Stating 
it another way, there appears to be only a 4 to 5% chance of a 
loss in a majority of the Cases. When the probability of 
success was reduced to 1 in 20 there was still a 62% chance of 
making a profit and, even when the number of wells drilled was 
reduced to 5 rather than 20, the probability of a profit was 
approximately 56% (the lowest value in any Case) . 

An inspection of Table 12-10 reveals that the substitu- 
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tion of the Beta distribution for discounted profit per barrel 
rather than the Triangular distribution has very little effect 
on the Best Estimate of Discounted Profit, E(P). The risk as 
measured by the Standard Deviation and the Coefficient of 
Variation is increased however, by about 4.4%. The minimum 
value of discounted profit with 90% Confidence Limits has been 
decreased by 4.34% and the maximum value has been increased by 
1.78%. Additionally, the range of values with 90% Confidence 
Limits has been increased with the Beta distribution by 2.24%. 

With regard to the Base data run at 10,000 iterations 
rather than 1000 iterations, it is seen from Table 12-10 that 
the value of E(P) only varies by 1.47% from the 1000 iteration 
value while the relative risk as measured by the Coefficient 
of Variation increases by 1.5% primarily due to the decrease 
in the absolute value of E(P). 

12.3 Effects of Parameter Changes on the Best Estimate of 
Discounted Profit and Coefficient of Variation 

The first question explored in this Section is the 
effect of a change in parameter value on the best estimate of 
discounted profit, E(P). This effect shall be assessed by 
calculating a Significance Ratio for each Case. The Signifi- 
cance Ratio is defined as the change in observed value (i.e. 
E(P)) to a change in the parameter value. In a conceptual 
sense the Significance Ratio measures the significance of a 
change in a parameter of the model on a value calculated by 
the model. We assume, of course, that all other parameters 
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(except the one being studied) remain the same. Mathematically, 
this can be expressed as: 



SR 



A Observed Value 
A Parameter 



( 1 ) 



With this Significance Ratio the relative importance of a 
change in various model parameters on the value being observed 
can be measured. For the profitability model it is possible, 
for example, to measure the relative importance of a change 
in dry hole costs, lowest profit per barrel, etc. on the best 
estimate of discounted profit for the total venture. For the 
analysis of the best estimate of discounted profit, Equation 
1 can be rewritten: 



SR 



A E (P) 

A Parameter 



( 2 ) 



The values to be inserted in Equation 2 are obtained from 
Tables 12-1A through 12-9A. As an example, from Table 12-1A 
the following Significance Ratio is obtained: 

For a change in parameter of -75.001 (Case 1) from the 
Base we have a commenserate reduction of E(P) of -76.77%. 
Inserting these values in Equation 2 we have: 



SR 



-76.77 

-75.00 



SR = 1.02 



For Case 2 we have the following relationship: 



SR 



-50.44 

-50.00 



SR = 1.0. 
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In a similar fashion the Significance Ratios can be calculated 
for each of the nine Cases involving a variation of the number 
of wells drilled as shown on Table 12-1 and 12-1A. If this is 
done Significance Ratios ranging from .81 to 1.18 with a mean 
value of 0.99 are obtained. 

Table 12-11 contains a summary of the Ranges and Mean 
Significance Ratios for each of the parameter sets contained on 
Tables 12-1 through 12-9 and Tables 12-1A through 12-9A. In 
addition, Table 12-11 contains a value of the Normalized 
Significance Ratio which reduces each Significance Ratio to a 
value on a 1-100 scale in proportion to its actual Significance 
Ratio. The Normalized Significance Ratio can be utilized to 
facilitate comparison of parameter sets 

TABLE 12-11 

SIGNIFICANCE RATIOS - E(P) 



Parameter Set 

n 

P 

Q at 5% 

Q at 95% 

Q (max) 

Q (min) 

C u 
e (a) 

e(b) 
e (c) 



Significance Ratios 

Range Mean Normalized Mean 



.81 - 


1.18 


.99 


.96 - 


1.10 


1.03 


.19 - 


.40 


.30 


.56 - 


.82 


.72 


.02 - 


.12 


.07 



0 0 

0 . 04 

0 .21 

0 .36 

0 .47 



96 

100 

29 

70 

07 

0 

04 

20 

35 



46 
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Several observations can be made from Table 12-11. 
First, it can be seen that a change in the number of wells 
drilled and probability of success are both almost equally 
significant as far as their effects on the value of E(P). In 
fact, there is virtually a 1 to 1 relationship between a change 
in the number of wells drilled and probability of success and a 
change in the best estimate of discounted profit. In other 
words, if the number of wells drilled is increased by 101, 
almost a 10% increase in discounted profit can be expected. 
Likewise, if the probability of success is decreased by 20% 
about a 20% decrease in the discounted profit can be expected. 

Also seen from Table 12-11 is that the next most 
significant parameter in the model is the estimated value of 
reserve size with 95% probability, Q at 95%. In fact, this 
parameter has a normalized significance of 70 which is approx- 
imately 2.5 times as critical as the lower estimate of the 
reserve size, Q at 5% probability, on the discounted profit. 

In a similar fashion it can be observed that variations in the 
values of the upper and lower boundary limits for reserve sizes 
have little or no significance on the discounted profit. 

Although one might intuitively feel that changes in dry 
hole costs would have a material effect on the discounted 
profit it is seen that it only has a normalized significance 
of 4 which, in essence, is relatively insignificant in compar- 
ison to the other parameters. 

In comparing the three parameters associated with the 
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discounted profit per barrel it is noted that, of the three, 
changes in the highest profit per barrel, e(c), are 2.3 times 
as significant as changes in the lowest value, e(a), and 1.3 
times as significant as changes in the most likely value, e(b), 
Also observed is that of the 10 parameters measured, the value 
of e(c) ranks fourth in significance. 



A similar analysis regarding risk can be made by 
computing the Significance Ratios for the Coefficient of 
Variation. In this case Equation 1 can be rewritten as: 

SR _ A_CV 

A parameter 

Table 12-12 contains the results of these calculations. 

TABLE 12-12 

SIGNIFICANCE RATIOS - CV 

Significance Ratios 



(3) 



Parameter Set 


Range 


Mean 


Normalized 


n 


.28 


- 1.36 


.67 


100 


P 


. 31 


- 1.55 


.64 


96 


Q at 51 


.04 


- .43 


.19 


28 


Q at 951 


. 01 


- .18 


.06 


9 


Q (max) 


. 06 


- .15 


.10 


15 


Q (min) 




0 


0 


0 


C 


.04 


- .05 


. 04 


6 


u 

e (a) 


.06 


- .11 


.09 


13 


e (b) 


0 


- .06 


. 03 


5 


e (c) 


.04 


- .08 


.06 


9 
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From Table 12-12 certain observations regarding the 
significance of a change in model parameters on the relative 
risk of the venture can be made. First, it is observed that a 
change in the number of wells drilled and in the probability 
of success both have about the same effect on the risk of the 
venture. It is interesting to note that the Normalized 
Significance Ratios for E(P) are exactly reversed for CV. In 
this regard, Table 12-11 indicated that a change in Q at 951 
had a greater significance than a change in Q at 51. Table 
12-12 reveals the opposite is true for the relative risk. 

Here a change in Q at 51 has a much greater significance than 
a change in Q at 95%. The same reversal is evident in the 
values of lowest and highest profit per barrel, e(a) and e(c). 

It must be remembered, of course, that the Signifi- 
cance Ratio for risk (coefficient of Variation) does not in- 
dicate whether the risk is increasing or decreasing with a 
parameter change, but only the relative amounts that it is 
changing. For the increase or decrease correlation we must 
go to Tables 12-1 through 12-10 and 12-1A through 12-9A. 

12.4 Convergence of the Expected Profit 

If the Monte Carlo simulation model is executed only 
once, only one value for the discounted profit is generated and 
there is no distribution of values that can be analyzed. Since 
simulation is being used because the model is too complex for 
other types of analyses, there is no basis for relying on the 
single result as a typical output of the model - there is no 
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way to tell whether the result is an extreme case or about the 
average . 



As the simulation process is repeated frequencies of 

occurrence are obtained from which the averages, the extremes, 

and what the shape of the distribution would be if the model 

had been executed infinitely many times can be inferred. In 

this regard, R.F. Barton [6] poses the following question: 

.... how reliable are the results from a limited 
number of executions? Two executions may give us a 
great deal more information than a single execution. 
One hundred iterations may give us hints of the extre- 
mes and shapes of the ultimate distributions. A thou- 
sand executions may perhaps give us results that app- 
roximate smooth curves. But a second thousand execu- 
tions may suggest smooth curves of different locations 
and shapes. A third thousand executions may generate 
results that look even different. On the other hand, 
the separate sets of one thousand executions may look 
so similar we feel no need for further executions. 



These possibilities raise a very practical question: 

How long to run a simulation model? One aspect of this quest- 
ion was analyzed by measuring the convergence of the discounted 
Expected Profit, E(P), averaged over the number of iterations 
run, with the discounted Expected Profit averaged over 1000 
iterations . 

As an example, if the value E(P) for the first iteration 
is $100.0 million and $50.0 million for the second iteration 
the discounted E(P) as determined by two iterations is $75.0 
million. We can then compare this value with the E(P) averaged 
over 1000 iterations, say $60.0 million, to determine the 
rapidity of convergence. If we were satisfied with an Expected 



128 



Profit that varied no more than ± 51 from the value which would 
be obtained with 1000 iterations (90% certainty) then by plot- 
ting or analyzing the data we could determine how many iter- 
ations were necessary to reach this point. It may turn out, in 
our example, that a value of $60.0 million ± $3.0 million was 
achieved after only 500 iterations. By performing this analysis 
with many different values of model variables, using the same 
random numbers, and each time determining the point beyond 
which the averaged E(P) varies by no more than ± 5% of the 1000 
iteration value, we may be able to infer the number of itera- 
tions necessary for 90% confidence limits for this particular 
model . 

This is, in fact, what has been done. For each set of 
model variables the iteration beyond which the value of E(P) 
does not vary by more than ± 5% from the value of discounted 
expected profit after 1000 iterations was measured. These 
values are contained on Table 12-13. Figures 12-9 through 
12-14 graphically depict this convergence process for various 
combinations of model parameters. 

In addition to this analysis, the model was run with 
the base data for a total of 10,000 continuous iterations to 
determine if any periodicities in the variables, as suggested 
by Barton [7], occur which would materially effect the value 
of the discounted Expected Profit after 10,000 iterations. 

Table 12-10 and Figure 12-15 contain the results of this test. 

It can be seen from Figure 12-15 that no periodicities occur 
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and the 10,000 iteration value of E(P) is, in fact, only 
slightly different from the 1000 iteration value. 

A visual inspection of the data contained on Table 
12-13 indicated that Expected Profit, E(P), does not converge 
with any regularity, either within 80% or 90% confidence limits. 
It is seen that even within a parameter set(i.e. number of 
wells drilled) the number of iterations vary greatly. The only 
significant fact which can be discerned from the data is that 
values of E(P) within 80% confidence limits can be achieved 
anywhere from 10 to 339 iterations and within 90% confidence 
limits anywhere from 63 to 710 iterations. If the highest 
and lowest iteration value is eliminated the range for 80% 
confidence limits is 18 to 242 iterations and the range for 
90% confidence limits is 68 to 596 iterations. 
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TABLE 12-13 

CONVERGENCE OF E(P) WITH NUMBER OF ITERATIONS 







WITHIN 


CONFIDENCE LIMITS 




Case 


Factor 

Being 


E(P) 

1000 


Within 80% Conf. 
Limits (± 10%) 

No. of 


Within 90% Conf. 
Limits (± 5%) 

No . of 


No. 


Varied 


Iterations 


I terations 


I terations 



Base 


- 


68.02 


117 


207 


1 


n 


15.80 


239 


502 


2 


n 


33.71 


242 


457 


3 


n 


39.49 


190 


440 


4 


n 


51.48 


186 


268 


5 


n 


60.02 


147 


195 


6 


n 


73.52 


65 


107 


7 


n 


83.12 


111 


261 


8 


n 


99.47 


136 


238 


9 


n 


115.37 


10 


63 


10 


p 


18.19 


339 


710 


11 


p 


33.75 


213 


268 


12 


p 


43.56 


204 


395 


13 


p 


81.08 


108 


255 


14 


p 


90.46 


41 


83 


15 


p 


100.50 


56 


167 


16 


p 


114.13 


42 


99 


17 


p 


125.12 


26 


68 


18 


P 


137.27 


40 


96 


19 


Q(5%) 


46.15 


120 


543 


20 


Q (5%) 


61.16 


117 


596 


21 


Q (5% ) 


63.24 


117 


206 


22 


Q (5% ) 


65.30 


117 


207 


23 


Q (5% ) 


69.34 


17 


129 


24 


Q (5% ) 


70.90 


17 


129 


25 


Q(5%) 


72.89 


17 


129 


26 


Q(5%) 


80.42 


42 


129 


27 


Q (95% ) 


23.45 


105 


129 


28 


Q (95% ) 


52.42 


17 


207 


29 


Q (95%) 


57.50 


17 


128 


30 


Q (9 5% ) 


62.64 


18 


207 


31 


Q (95%) 


72.63 


117 


129 


32 


Q (95%) 


77.14 


117 


207 


33 


Q (9 5% ) 


81.69 


118 


207 


34 


Q (9 5% ) 


97.46 


119 


473 



Casi 

No 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 
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TABLE 


12-13 Continued 




Factor 

Being 

Varied 


E(P) 

1000 

I terations 


Within 801 Conf. 
Limits (± 101) 

No. of 
Iterations 


Within 90% Conf. 
Limits (± 5%) 

No. of 
Iterations 



Q max 


64.35 


43 


129 


Q max 


69.29 


18 


305 


Q min 


68.02 


117 


207 


Q min 


68.02 


117 


207 


Q min 


68.18 


22 


129 


Cu 


69.30 


18 


129 


Cu 


68.87 


18 


129 


Cu 


67.17 


117 


207 


Cu 


65.47 


118 


207 


Cu 


62.07 


117 


265 


Cu 


53.58 


119 


539 


e (a) 


60.99 


117 


207 


e (a) 


63.80 


117 


207 


e (a) 


65.21 


117 


129 


e (a) 


66.61 


117 


207 


e (a) 


69.43 


18 


207 


e (a) 


70.83 


19 


129 


e (a) 


72.24 


18 


129 


e (b) 


60.74 


117 


129 


e (b) 


62.97 


22 


129 


e (b) 


65.50 


22 


129 


e (b) 


70.54 


117 


207 


e (b) 


73.07 


117 


207 


e (b) 


75.31 


117 


126 


e (c) 


60.54 


117 


207 


e (c) 


61.65 


117 


207 


e (c) 


64.70 


117 


207 


e (c) 


71.35 


18 


129 


e (c) 


74.39 


22 


129 


e (c) 


77.71 


22 


129 


e (c) 


99.87 


22 


129 


Base 


67.98 


166 


239 
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Figure 12-11 Convergence of E(P) 
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Figure 12-12 Convergence of E(P) - Case 45; C = $1.0 Million 
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Figure 12-13 Convergence o£ ECP) " Case 52; e(a) - $.65 
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Figure 12-14 Convergence of E(P) - Beta Distribution 
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Figure 12-15 Convergence of E(P) - Base, 10,000 Iterations 
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12.5 Changes in Return and Risk with the Number of Wells 
Drilled 

Since the risk of the venture has been measured by the 
coefficient of the variation it may be of interest to ascertain 
whether there is a mathematical relationship between the risk 
and the number of wells drilled. To determine whether a mathe- 
matical relationship exists we start with Equation 7, Chapter V 



CV 





Now, if p is held constant it can be seen that the Coefficient 
of Variation should be proportional to the number of wells 
drilled as follows: 



cv ~7Z 


(4) 


CV x /n = Constant 


(5) 



Table 12-14 contains a test of this proportionality hypothesis 
using the data contained on Table 12-1 



TABLE 12-14 

PROPORTIONALITY TEST - CV to n 
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n 


/n 


Coeff. of Variation 


Constant 


5 


2 . 23 


1.59 


3.54 


10 


3.16 


1.11 


3.52 


12 


3.46 


1.07 


3.70 


15 


3.87 


0.92 


3.56 


18 


4.24 


0.84 


3.56 


20 


4.47 


0.79 


3.53 


22 


4.69 


0.75 


3.52 


25 


5.00 


0.71 


3.55 


30 


5.48 


0.64 


3.50 


35 


5.92 


0.62 


3.67 



It can be seen from Table 12-14 that the Coefficient of Vari- 
ation is indeed proportional to l//n. Thus, if the number of 
wells is increased by 4 times the risk will be reduced by a 
factor of 2. 

It is also interesting to note that the Best Estimate 
of Discounted Profit per well, or E(P)/n, remains fairly con- 
stant while the risk, as measured by the Coefficient of Vari- 
ation, decreases with l//n. This can be seen on Table 12-15. 
This would tend to indicate that while the best estimate of 
return per well remains fairly constant with the number of 
wells drilled the risk associated with that return deminishes 
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as we increase the number drilled. 



TABLE 12-15 

BEST ESTIMATE OF PROFIT PER WELL AND RISK 



Case 


Number of 
Wells Drilled 
(n) 


Best Estimate 
of Discounted 
Profit- E (P) 


E (P) 
n 


Coeff. of 
Variation 
CV 


Constant 

1 

/n 


1 


5 


15.80 


3.16 


1.59 


3.54 


2 


10 


33.71 


3.37 


1.11 


3.52 


3 


12 


39.49 


3.29 


1.07 


3.70 


4 


15 


51.48 


3.43 


0.92 


3.56 


5 


18 


60.02 


3.33 


0.84 


3.56 


Base 


20 


68.02 


3.40 


0.79 


3.53 


6 


22 


73.52 


3.34 


0.75 


3.52 


7 


25 


83.12 


3.32 


0.71 


3.55 


8 


30 


99.47 


3.32 


0.64 


3.50 


9 


35 


115.37 


3.29 


0.62 


3.67 



CHAPTER XIII 



CONCLUSIONS 



13.1 General 

The following general conclusions are drawn from the 
results of this study: 

1. The model developed herein provides a suitable 
vehicle for measuring the sensitivity of the parameters affect- 
ing the profitability of an exploratory oil drilling venture. 

2. The parameters in the profitability model should 
include not only single best estimate values but, where app- 
ropriate, also a probability distribution of values. 

3. The probability distribution best representing the 
number of successes is the Binomial distribution, especially 
when the probability of success and number of wells drilled are 
in the order of magnitude normally encountered in an exploratory 
drilling program. For probabilities of success approaching 50% 
and large numbers of wells drilled, the Normal or Gaussian dis- 
tribution can closely approximate the Binomial distribution. 

4. The reserve size discovered by a successful well is 
best represented by a Lognormal distribution. The values can 
be calculated manually using points scaled from a straight line 
drawn on log-probability paper between values of reserve size 
which can be expected with 5 and 95% probability. The lognormal 
distribution can be simulated in computer applications with 
excellent correlation using the Lognormally Distributed Variable 
Generator contained in this study. 
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5. The profit per barrel received over a period of 
years from oil produced must be discounted to reflect the time 
value of money. 

6. The Triangular distribution function adequately 
serves to represent the discounted profit per barrel. The Beta 
distribution function provides only slightly different values 
of discounted profit per barrel when tested using the same 
input data used in the Triangular distribution. Both the 
Triangular and Beta distributions can be simulated in computer 
applications using the Generators contained herein. 

7. The complexity of the distributions in the model 
precludes manual computation of profitability and sensitivity 
studies thereby necessitating the use of the Monte Carlo 
technique and a digital computer. 

13.2 Sensitivity Studies 

In the sensitivity studies, where individual parameters 
in the model were varied, with the remaining parameters being 
held constant, the following conclusions were drawn: 

1. The one single observed value which fluctuated the 
greatest (on a percentage basis) was the minimum value of 
discounted profit within 90% Confidence Limits. 

2. Except in isolated cases, the chance if a profit 
greater than zero in all tests remained fairly constant at 95 
to 96%. Conversely, there remained a fairly constant chance 
of 4 to 5% of incurring a loss. 
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3. A subsitiution of the Beta distribution for dis- 
counted profit per barrel rather than the Triangular distri- 
bution makes very little effect on the best estimate of dis- 
counted profit for the venture or the risk as measured by the 
coefficient of variation. 

4. There appears to be no periodicities in the para- 
meters of the model and the best estimate of discounted profit 
calculated after 10,000 executions varies only slightly from 
the value calculated after 1,000 iterations. 

5. It is possible to calculate values of the best 
estimate of discounted profit for the venture within 80% 
confidence limits (± 10% of the 1,000 iteration value of E(P)) 
by executing the model between 18 and 242 times; and within 
90% confidence limits by executing the model between 68 and 
596 iterations. Further, there appears to be no regular 
convergence of E(P) with the number of iterations performed. 

6. A Significance Ratio, relating a change in observed 
value calculated by the model to a change in parameter value, 
provides an excellent tool for comparing the effects of changes 
of different model parameters, especially when the Significance 
Ratio is normalized on a 1 to 100 scale. 

7. Using the Significance Ratios to relate parameter 
changes to discounted venture profitability, E(P) it is con- 
cluded that with the parameter values used: 

a. Changes in probability of success and number 
of wells drilled have almost equal significance on changes in 



venture profitability. 
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b. Changes in the lower range of reserve size have 
only about 30% of the effect of changes in p or n on E(P). 
However, changes in the upper range have a 70% significance. 

The values selected at lower and upper boundaries for reserve 
sizes in the model have little or no effect on the value or E(P). 

c. Changes in the highest value of discounted 
profit per barrel have more than twice the significance as 
changes in the lowest values of discounted profit per barrel; 
but even then, have only 46% of the significance as changes 
in n or p on E (P) . 

d. Changes in dry hole cost have relatively 
little effect on E(P). 

8. Using the Significance Ratios to relate parameter 
changes to the relative risk of the venture as measured by the 
coefficient of variation, CV , it is concluded that with the 
parameter values used: 

a. Changes in the number of wells drilled, n, and 
the probability of success, p, have about the same effect on 
the risk of the venture. 

b. Changes in the lower range of reserve size are 
next in significance. Changes in this parameter are more than 
3 times as significant as changes in the upper ranges of 
reserve sizes. 

c. The value selected as the upper boundary for 
reserve size will effect the risk in the venture but the lower 
boundary limit will not. 



d. Changes in dry hole costs do little to change 
the risk of the venture. 

e. Changes in the lowest value of discounted 
profit per barrel is about 1.5 times as significant as changes 
in the highest value of discounted profit per barrel but even 
then, has only a significance of 13% as compared with the risk 
associated with changes in the number of wells drilled. 

9. While the return per well remains fairly constant 
as the number of wells increase, the relative risk of the 
venture decreases as the square root of the number of wells 
drilled. 



The seemingly large number of findings drawn from the 
results of this study are indicative of the value of this type 
of analysis. It is recognized, of course, that not all the 
conclusions listed herein will be considered of equal import- 
ance by the reader - it would be rare if they did. The fact 
that these findings could be drawn at all, however, leads to 
the overall conclusion that sensitivity analysis of this type 
provides a valuable tool for the decision-maker. 



CHAPTER XIV 



RECOMMENDATIONS 

It is in the light of the problems facing most decision- 
makers, namely the "What if...?" questions cited in Chapter XI 
and the contention so lucidly expressed by William T. Morris 
[35] that at some point we must stop gathering data and re- 
fining our information and start recommending action, that the 
profitability model is ideally suited. 

It is therefore recommended that a model of this type 
be used not to obtain absolute data on whether a profit will 
be made for the venture (or some other observation of this 
type) since other methods of analysis may be better suited for 
this purpose. Instead, it is recommended that the model be 
used to assess the relative sensitivity of the results to para- 
meter changes in the model. With this knowledge of parameter 
sensitivity we can now devote time and effort to refining only 
those parameters which have a significant affect on the outcome. 

For example, if we were concerned with the best estimate 
of discounted profit for the venture we would probably devote 
our efforts to refining our estimates of the probability of 
success, the upper range of reserve size and the highest value 
of discounted profit per barrel (assuming, of course, that the 
number of wells drilled is fixed). We would not, in this 
example, concern ourselves with such items as the dry hole 
costs, the lower range of reserve sizes and the like. 
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In this way, using the sensitivity analysis approach, 
the decision-maker can achieve the maximum return for the 
effort he can afford to expend. 



APPENDIX A 
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THE LOGNORMALLY DISTRIBUTED VARIABLE GENERATOR 

The FORTRAN IV Computer program for lognormally 
distributed deviates requires two different, uniformally 
distributed random fractions to generate one normal variable. 
With two uniformally distributed fractions, RA and RB, one 
generates a variable, V, which is normally distributed about 
a mean of 0, with a standard deviation of 1, thus, 

V = (-2 In RA)’ 5 x COS 2tt (RB) 

To convert V to a variable QNORM, which is normally distrib- 
uted about a mean of TMU, with a standard deviation SD, we 
simply use the relationship, 

QNORM = (V) SD + TMU 

By multiplying the natural logarithm base, e, to the normal 
deviate power, a lognormally distributed deviate, QLOGN, is 
generated . 

QLOGN = e Q N0RM 

Once the Lognormally distributed deviate is generated 
it is tested against the lower and upper limils established as 
inputs to the program. If the deviate is less than the minimum 
value of QMIN, it is set to read QMIN. If it is greater than 
the maximum value of QMAX, it is set to read QMAX. 

It will be noted that the program requires a value of 
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mean and standard deviation. These values can be obtained 
from the median value of Q at W( ^Q) = 50%, and the value of 
Q at W( ^Q) = 84% (or, y m + Oy) as follows: 



From Equation 3, 

y m (median) = y Q (mean) 

From Equation 4B, 

ln x m = y 0 

Qmedian = x m 
and , TMU = y Q 

then , 

TMU = log e (Qmedian) 
or 

TMU = ALOG (QMEDIAN) 



y 0 = y m 



y m + a y 
y m + a y 

or, 

y m + a y 

and , 

m x m 

Therefore , 

(>m + °y) - ym 

or, 



y G + °y 

ln Q at W(<Q) = 84% 



ln Q (84) 



y m 

ln Q (84) - In (Q m ) 



SD = ALOG (Q84) 



ALOG (QMEDIAN) 
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The FORTRAN IV Computer program incorporating these concepts 
is listed as follows: 

PROGRAM QLOGN (INPUT , OUTPUT) 

READ 100, QMEDIAN , Q84, QMAX, QMIN 
SD = ALOG(Q84) - ALOG (QMEDIAN) 

TMU = ALOG (QMEDIAN) 

RA = RANF (0) 

RB = RANF (0) 

V = (- 2 . 0*ALOG (RA) ** . 5* COS (2 .0*3.1416*RB) 
QNORM = V*SD + TMU 
QLOGN = EXP (QNORM) 

IF (QLOGN - QMIN) 6, 7, 8 

6 QLOGN = QMIN 
GO TO 7 

8 IF (QLOGN - QMAX) 7, 7, 9 

9 QLOGN = QMAX 

7 CONTINUE 

PRINT 120, QLOGN 
100 FORMAT (10X, 4F10 . 3) 

120 FORMAT (20X,F10.3) 



END 



APPENDIX B 
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THE TRIANGULAR DISTRIBUTION VARIABLE GENERATOR 



The FORTRAN IV Computer program for deviates with a 
Triangular distribution requires one uniformally distributed 
random function, RANF(O). This uniformally distributed random 
function will assume either the value of W(x < x^) , - WXLTXI - 
or W(x > x^) - WXGTXI - depending on whether the calculated 
value of x^(XI) is equal to or less than the inputed value of b. 

The FORTRAN IV program incorporating these concepts is 



listed as follows: 

PROGRAM TRIA (INPUT , OUTPUT) 

THIS PROGRAM CALCULATES THE BEST ESTIMATE OF PROFIT/BBL 
ASSUMING A * TRIANGULAR * DISTRIBUTION OF PROBABILITY. 

A = MINIMUM PROFIT/BBL EXPECTED 
B = MOST LIKELY VALUE OF PROFIT/BBL EXPECTED 
C = MAXIMUM PROFIT/BBL EXPECTED 

READ 112, A, B, C, 

CALL RANF(O) 

WXLTXI = RANF(O) 

WXGTXI - 1.0 - WXLTXI 

XI = A + SQRT (WXLTXI* (B-A) * (C-A) ) 

IF (XI.LE.B) 55,50 

50 XI = C - SQRT (WXGTXI* (C-B)*C-A)) 

55 PROFIT = XI/100. 

PRINT 300, PROFIT 
112 FORMAT (10X , 3F10 . 2) 

300 FORMAT (15X,F4.2) 

END 



APPENDIX C 



156 



THE BETA DISTRIBUTION VARIABLE GENERATOR 



The FORTRAN IV Computer program for deviates with a 



Beta distribution is listed as follows: 



PROGRAM BETA ( INPUT , OUTPUT) 

THIS PROGRAM CALCULATES THE BEST ESTIMATE OF PROFIT/BBL 
ASSUMING A * BETA * DISTRIBUTION OF PROBABILITY. 

A = MINIMUM PROFIT/BBL EXPECTED 
B = MOST LIKELY VALUE OF PROFIT/BBL EXPECTED 
C = MAXIMUM PROFIT/BBL EXPECTED 

READ 100, A, B, C 
XMU = (4.*B+A+C)/6. 0 
XVAR = (C -A) * * 2/ 36 . 0 
BMEAN = (XMU- A)/ (C-A) 

BVAR = XVAR/ (C-A) ** 2 

XK1 = BMEAN* (BMEAN* (1. 0 - BMEAN) /BVAR- 1 . 0) 

XK2 = XK1* ( (1 . 0- BMEAN) /BMEAN) 

CALL GAMMARN (XK1,GAM1) 

CALL GAMMARN (XK2,GAM2) 

BETA = (GAM1/ (GAM1 + GAM2))*(C-A) +A 
PROFIT = BETA/100 . 

PRINT 300, PROFIT 
100 FORMAT (10X, 3F10 . 2) 

300 FORMAT (15X,F4.2) 

END 

SUBROUTINE GAMMARN (TK , GAM) 

GAMMA =1.0 
K1 = TK 
TK1 = K1 
Y1 = RANF(O) 

IF (Yl- (TK-TK1) ) 10, 10, 20 
10 K1 = K1 + 1 
20 DO 30 I = 1 , K1 
Y2 = RANF(0) 

30 GAMMA = GAMMA* Y 2 
GAM = -ALOG (GAMMA) 

RETURN 

END 
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COMPUTER PRINTOUT FOR BASE DATA 



PROFIT DISTRIBUTION OF AN EXPLORATORY DRILLING PROGRAM 



NO. 


OF WELLS 


prob of 


DRY HOLE 






DRILLED 




SUCCESS 


COST 


gmedian 


Q84 QMAX 


QMIN 




20 




.150 


.150 


20.00 


46.00 140.00 


1.00 




THE 


TRIANGULAR 


DISTRIBUTION FOR 


DETERMINATION OF 


THE 




PROFIT PER BARREL HAS 


THE FOLLOWING EXTREMES IN 


CENTS - 




A 




50.00 


B = 


88.00 


C = 115.00 




ITERATION 


NUMBER OF 


FIELD 


PROFIT PER 


ITERATION BEST ESTIMATE 






SUCCESSES 


SIZES 


BARREL 


PROFIT 


OF PROFIT 


1 




- 


2 


52.90 


.83 


41.17 


41.17 


2 






3 


64.03 


1.01 


61.90 


51.53 


3 






0 


0.00 


0.00 


-3.00 


33.36 


A 






4 


114.52 


l.o2 


114.13 


53.55 


5 






5 


189.57 


.79 


148.40 


72.52 


6 






0 


0.00 


0.00 


-3.00 


59.93 


7 






3 


34.12 


.87 


27.05 


55.24 


8 






5 


254.88 


1.01 


254.59 


80.15 


9 






3 


198.18 


.81 


157.93 


88.80 


10 






2 


32.61 


1.01 


30.36 


82.95 


11 






0 


0.00 


0.00 


-3.00 


75.14 


12 






2 


38.72 


.93 


33.17 


71.64 


13 






5 


179.08 


.93 


164.57 


78.79 


14 






3 


95.99 


.57 


51.88 


76.87 


15 






3 


166.86 


.61 


98.92 


78.34 


16 






3 


129.04 


.65 


81.67 


78.55 


17 






2 


28.20 


.86 


21.46 


75.19 


18 






3 


54.65 


.80 


41.00 


73.29 


19 






3 


73.73 


.95 


67.13 


72.97 


20 






1 


85.62 


.97 


80.58 


73.35 


21 






5 


127.41 


.81 


100.96 


74.66 


22 






1 


3.84 


.92 


.69 


71.30 


23 






2 


40.45 


.79 


29.10 


69.46 


24 






4 


59.44 


.87 


49.05 


68.61 


25 






2 


53.91 


.81 


41.17 


67.52 


26 






3 


49.27 


.66 


30.22 


66.08 


27 






2 


26.85 


.70 


16.03 


64.23 


28 






5 


96.30 


.73 


67.94 


64.36 


29 






2 


27.60 


.74 


17.80 


62.75 


30 






4 


187.40 


.99 


183.55 


66.78 


31 






2 


137.95 


.64 


85.37 


67.38 


32 






4 


83.15 


.82 


65.70 


67.33 


33 






2 


107.08 


.82 


84.63 


67.85 


34 






0 


0.00 


0.00 


-3.00 


65.77 


35 






4 


84.91 


.85 


69.90 


65.89 


36 






2 


38.40 


.94 


33.52 


64.99 


37 






2 


25.44 


1.02 


23.35 


63.86 


38 






4 


125.78 


.73 


89.62 


64.54 


39 






1 


6.76 


.78 


2.41 


62.95 


40 






5 


110.73 


.89 


95.83 


63.77 


41 






0 


0.00 


0.00 


-3.00 


62.14 


42 






4 


116.96 


1.09 


124.90 


63.63 



43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 
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6 229.75 .91 

1 32.97 .97 

25.19 1.07 

0.00 0.00 

181.65 1.02 

84.76 .87 

138.59 .59 

160.84 1.04 

91.12 1.04 

50.01 .71 

30.00 .67 

43.41 .55 

129.75 .94 

19.93 .79 

63.36 .92 

72.45 .87 

110.63 .63 

100.93 .66 

194.06 .94 

170.74 .70 

0.00 0.00 

74.02 .86 

39.45 .87 

89.92 .96 

112.67 .84 

172.22 .88 

51.86 .99 

47.19 .84 

63.76 .91 

76.24 .58 

23.86 .80 

110.39 .80 

20.76 .84 

95.34 .93 

38.74 .87 

23.97 .76 

32.79 .79 

50.87 .89 

91.19 .73 

96.51 .63 

119.23 .91 

0.00 0.00 

184.12 .73 

31.18 .78 

51.86 .86 

80.81 .58 

109.49 .95 

55.49 .57 

19.07 1.12 

29.24 .94 

88.87 .70 

23.84 .88 

83.49 .67 

88.52 1.06 

156.23 .84 

0.00 0.00 

50.68 1.05 

42.56 .93 

0.00 0.00 

55.47 .87 

88.88 .91 

34.29 .89 



206.55 


66.96 


29.02 


66.10 


24.20 


65.17 


- 3.00 


63.68 


183.29 


66.23 


70.92 


66.33 


79.57 


66.60 


165.66 


68.58 


92.52 


69.05 


32.64 


68.35 


17.45 


67.39 


20.90 


66.53 


119.35 


67.49 


12.84 


66.51 


55.84 


66.32 


60.17 


66.22 


66.80 


66.23 


64.42 


66.20 


181.04 


68.08 


117.80 


68.88 


- 3.00 


67.74 


60.82 


67.63 


31.43 


67.08 


83.53 


67.32 


91.91 


67.69 


150.10 


68.90 


48.86 


68.61 


36.78 


68.16 


55.78 


67.98 


41.86 


67.62 


16.45 


66.92 


85.40 


67.17 


14.53 


66.47 


85.91 


66.72 


30.88 


66.26 


15.85 


65.61 


23.31 


65.08 


42.79 


64.80 


64.26 


64.79 


58.39 


64.71 


106.21 


65.21 


- 3.00 


64.40 


132.41 


65.20 


21.46 


64.69 


41.95 


64.43 


44.18 


64.20 


101.59 


64.62 


28.98 


64.23 


18.70 


63.72 


24.86 


63.30 


59.78 


63.26 


18.05 


62.78 


53.13 


62.68 


92.78 


63.00 


129.02 


63.68 


- 3.00 


63.00 


50.79 


62.87 


37.08 


62.61 


- 3.00 


61.96 


45.67 


61.81 


78.40 


61.97 


27.59 


61.64 



105 

106 

107 

108 

109 

no 

in 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

133 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 






124.19 


.68 


82.68 


61.84 


82.08 


1.13 


90.48 


62.11 


92.60 


.65 


57.52 


62.06 


141 .82 


.75 


103.75 


62.45 


101.35 


.87 


85.45 


62.66 


45.26 


.77 


32.31 


62.38 


76.50 


.77 


56.34 


62.33 


27.12 


.99 


24.04 


61.99 


18.92 


.96 


15.36 


61.58 


83.56 


. 83 


67.61 


61.63 


31 .66 


.84 


24.00 


61.30 


85.46 


.64 


52.03 


61 .22 


166.72 


.84 


137.93 


61.88 


199.84 


.68 


133.63 


62.49 


118.33 


.94 


108.42 


62.87 


150.97 


.85 


125.42 


63.39 


0.00 


0.00 


- 3.00 


62.84 


256.55 


.60 


151.98 


63.57 


199.11 


.93 


184.32 


64.56 


111.75 


.84 


91.60 


64.77 


25.56 


.80 


17.70 


64.40 


0.00 


0.00 


- 3.00 


63.86 


57.89 


.73 


39.42 


63.67 


83.49 


1.03 


83.96 


63.83 


245.33 


.95 


229.78 


65.11 


194.04 


.95 


181.68 


66.01 


108.37 


1.02 


108.45 


66.34 


72.57 


.54 


36.64 


66.11 


18.69 


1 . M 


16.22 


65.74 


59.56 


.86 


48.53 


65.61 


200.74 


.92 


182.97 


66.48 


128.45 


.95 


119.53 


66.87 


3.87 


.88 


.55 


66.38 


231.08 


.74 


169.64 


67.13 


46.63 


.67 


28.36 


66.85 


31.68 


.66 


18.03 


66.50 


52.13 


.88 


43.28 


66.34 


100.82 


1.06 


104.06 


66.60 


27.43 


1.10 


27.35 


66.33 


0.00 


0 c 00 


- 3.00 


65.85 


8.60 


.79 


3.96 


65.42 


246.11 


.93 


227.67 


66.53 


39.53 


.71 


25.60 


66.25 


209.61 


.99 


205.16 


67.19 


84.89 


.66 


53.98 


67.10 


33.16 


.75 


22.37 


66.81 


37.88 


.85 


29.71 


66.56 


5.58 


.65 


.80 


66.13 


266.54 


.71 


187.86 


66.92 


47.87 


.74 


32.89 


66.70 


30.16 


.96 


26.68 


66.44 


196.48 


.79 


152.41 


67.00 


161.99 


1.00 


159.09 


67.58 


89.22 


.83 


71.30 


67.61 


40.85 


.79 


29.58 


67.37 


135.26 


.80 


105.98 


67.61 


31.66 


1.06 


30.89 


67.38 


95.17 


.91 


84.64 


67.49 


79.70 


.84 


64.64 


67.47 


5.40 


.99 


2.50 


67.07 


90.85 


.85 


74.85 


67.12 


66.02 


.73 


45.18 


66.99 



167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 



3 


84.28 


1.05 


86.29 


67.10 


3 


87.96 


.89 


75.31 


67.15 


4 


48.41 


.96 


44.18 


67.02 


4 


47.52 


.95 


42.63 


66.87 


3 


49.62 


• 80 


37.04 


66.70 


4 


99.74 


.92 


89.48 


66.83 


0 


0.00 


0.00 


- 3.00 


66.43 


4 


57.47 


.99 


54.24 


66.36 


4 


163.81 


.75 


120.32 


66.66 


3 


155.76 


1.00 


153.48 


67.16 


4 


185.30 


• 68 


124.21 


67.48 


0 


0.00 


0.00 


- 3.00 


67.08 


0 


0.00 


0.00 


- 3.00 


66.69 


5 


108.78 


.97 


102.99 


66.89 


4 


82.87 


.78 


61 .96 


66.87 


2 


68.33 


• 81 


52.39 


66.79 


2 


43.93 


.99 


40.75 


66.65 


2 


92.28 


.62 


54.36 


66.58 


2 


49.77 


.72 


32.94 


66.40 


1 


10.66 


.82 


5.92 


66.07 


4 


101.08 


• 80 


78.67 


66.14 


3 


84.32 


1.00 


81.58 


66.22 


3 


116.27 


.94 


106.32 


66.43 


3 


46.74 


.77 


33.63 


66.26 


4 


63.46 


.78 


47.11 


66.16 


1 


69.48 


.75 


49.50 


66.07 


2 


33.12 


• 86 


25.72 


65.87 


0 


0.00 


0.00 


- 3.00 


65.51 


3 


127.19 


.85 


105.54 


65.72 


4 


146.63 


.94 


135.76 


66.07 


2 


48.18 


.70 


31.24 


65.90 


4 


81.22 


1.03 


81.18 


65.97 


4 


76.37 


.78 


57.05 


65.93 


2 


46.29 


.76 


32.42 


65.76 


1 


49.99 


.54 


24.21 


65.55 


2 


56.42 


.70 


36.53 


65.41 


1 


10.27 


• 81 


5.52 


65.12 


2 


33.89 


.89 


27.39 


64.93 


1 


12.06 


.78 


6.59 


64.65 


3 


67.32 


.87 


55.92 


64.60 


4 


146.58 


.90 


128.85 


64.91 


3 


78.10 


.70 


51.99 


64.85 


8 


221.96 


.71 


155.14 


65.28 


3 


32.79 


.78 


23.02 


65.08 


7 


179.39 


.62 


109.95 


65.30 


4 


102.93 


.65 


65.00 


65.29 


5 


182.99 


.59 


106.45 


65.49 


2 


23.19 


.90 


18.16 


65.27 


4 


179.31 


.85 


149.58 


65.66 


3 


40.05 


.89 


33.09 


65.51 


4 


68.35 


.79 


51.88 


65.44 


4 


94.78 


.90 


82.79 


65.52 


5 


171.29 


.97 


164.52 


65.98 


0 


0.00 


0.00 


- 3.00 


65.66 


5 


110.02 


.92 


98.94 


65.81 


2 


42.10 


.68 


26.12 


65.63 


5 


142.72 


.95 


132.74 


65.94 


5 


132.25 


.82 


105.64 


66.11 


2 


60.06 


.70 


39.54 


65.99 


2 


16.24 


.70 


8.69 


65.74 


3 


65.40 


.94 


58.74 


65.71 


3 


72.30 


1.04 


72.30 


65.74 



229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 

281 

282 

283 

284 

285 

286 

287 

288 

289 

290 



►-tv mow uuii\>ro.C'U>~j.C'-C'*-r\jr\jr\Ju>uju)o^iui*-i- ui-t'*-uro^uirur\> H-uroM-t'row^uwrorooui^ Hj-uroMO'u^ -t- 



72.06 


.78 


53.48 


65.69 


128.08 


.89 


111.51 


65.89 


45.47 


.67 


27.83 


65.72 


200.29 


1.01 


200.10 


66.30 


38.19 


.75 


25.98 


66.13 


67.98 


• 86 


56.05 


66.08 


43.40 


.76 


30.58 


65.93 


121.19 


.84 


99.09 


66.07 


30.96 


.73 


19.73 


65.88 


70.22 


.61 


40.32 


65.77 


90.32 


.66 


57.49 


65.74 


0.00 


0.00 


- 3.00 


65.45 


84.62 


.92 


75.16 


65.49 


63.96 


.82 


49.89 


65.42 


119.11 


.63 


72.76 


65.45 


75.18 


.90 


64.79 


65.45 


163.88 


.85 


136.17 


65.74 


62.32 


.86 


51.07 


65.68 


43.48 


• 80 


32.27 


65.55 


82.25 


.96 


76.80 


65.59 


51.07 


.74 


35.23 


65.47 


35.34 


.59 


18.11 


65.28 


62.71 


.98 


59.10 


65.26 


15.34 


.96 


11.95 


65.04 


113.74 


.80 


80.76 


65.14 


62.76 


.62 


36.00 


65.02 


183.22 


.76 


137.64 


65.31 


51.97 


1.08 


53.87 


65.26 


32.48 


.70 


19.95 


65.09 


95.57 


1.01 


93.89 


65.20 


24.58 


. 86 


18.21 


65.02 


100.22 


.87 


85.22 


65.09 


126.84 


.96 


119.65 


65.30 


21.42 


.95 


17.54 


65.12 


14.90 


.85 


9.85 


64.91 


81.39 


.90 


71.00 


64.93 


192.19 


1.04 


197.50 


65.43 


0.00 


0.00 


- 3.00 


65.18 


166.32 


.75 


121.45 


65.39 


115.73 


.76 


85.21 


65.46 


237.11 


.98 


230.96 


66.08 


82.53 


.93 


73.67 


66.11 


27.87 


• 86 


21.19 


65.94 


28.30 


.77 


19.11 


65.77 


10.30 


.77 


5.10 


65.55 


136.50 


.84 


111.76 


65.71 


97.33 


.86 


80.92 


65.77 


110.15 


.77 


83.08 


65.83 


117.57 


.73 


83.53 


65.90 


56 • 03 


.70 


36.79 


65.79 


66.87 


.75 


47.18 


65.72 


104.59 


.65 


65.00 


65.72 


81.22 


.74 


57.78 


65.69 


56.46 


1.10 


59.64 


65.67 


102.32 


1.10 


109.98 


65.83 


0.00 


0.00 


- 3.00 


65.59 


135.46 


.80 


105.93 


65.73 


150.39 


.57 


84.00 


65.79 


20.62 


.65 


10.58 


65.60 


170.68 


1.05 


176.30 


65.98 


25.73 


.69 


15.04 


65.81 


39.28 


1.05 


38.54 


65.71 



291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 
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54.14 


67.11 


105.29 


.76 


77.13 


67.13 


29.82 


.80 


21.06 


67.06 


40.81 


.95 


36.16 


67.01 


255.32 


.88 


221.15 


67.25 


65.67 


.86 


53.67 


67.23 


33.75 


.61 


17.90 


67.15 


35.03 


1.01 


32.47 


67.10 


216.11 


.64 


135.94 


67.20 


129.31 


.69 


86.98 


67.23 



663 

664 

665 

666 

667 

668 

669 

670 

671 

672 

673 

674 

675 

676 

677 

678 

679 

680 

681 

682 

683 

684 

685 

686 

687 

688 

689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

706 

707 

708 

709 

710 

711 

712 

713 

714 

715 

716 

717 

718 

719 

720 

721 

722 

723 

724 
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4 106.15 1.03 

2 32.78 .84 

7 152.33 1.03 

4 83.78 .92 

3 134.85 .97 

2 53.67 .80 

2 78.43 .86 

4 55.03 1.00 

3 114.17 .80 

2 40.13 .53 

4 69.77 .89 

1 27.81 .82 

2 19.20 .88 

3 56.07 .97 

4 135.30 .82 

3 66.58 .68 

4 50.25 .85 

5 145.83 .84 

1 56.21 .83 

131.39 .87 

0.00 0.00 

152.17 .73 

181.01 .82 

96.48 1.00 

56.78 .88 

145.81 .72 

110.81 .72 

36.83 1.00 

180.77 .88 

1 28.27 .88 

4 88.48 .92 

3 97.26 .72 

2 32.64 .84 

4 146.85 .71 

2 117.85 .90 

4 93.86 .76 

0 0.00 0.00 

2 51.79 .75 

5 78.39 .64 

3 84.02 .77 

4 87.68 .95 

3 97.17 .55 

2 26.85 .78 

2 42.73 .87 

2 69.03 .81 

1 30.12 .89 

3 37.92 .76 

3 75.44 1.01 

4 79.03 1.01 

4 92.93 .86 

3 114.75 1.11 

2 34.27 .84 

5 71.37 .77 

4 160.63 .69 

7 156.55 .59 

1 42.96 .97 

7 186.02 .84 

4 71.17 .54 

4 98.80 1.05 

2 16.43 .73 

4 144 .47 1.11 

1 73.36 .61 



106.92 


67.29 


24.79 


67.23 


154.63 


67.36 


74.91 


67.37 


128.37 


67 .46 


40.09 


67.42 


64.49 


67.42 


52.62 


67.39 


88.82 


67.43 


18.45 


67.35 


59.79 


67.34 


20.07 


67.27 


14.20 


67.19 


51.60 


67.17 


108.31 


67.23 


42.48 


67.19 


40.38 


67.16 


120.27 


67.23 


43.59 


67.20 


111.34 


67.26 


- 3.00 


67.16 


108.42 


67.22 


145.26 


67.33 


94.43 


67.37 


47.34 


67.35 


102.90 


67.40 


77.17 


67.41 


34.23 


67.36 


157.78 


67.49 


22.08 


67.43 


78.93 


67.44 


67.83 


67.45 


24.56 


67.38 


101.38 


67.43 


103.15 


67.48 


68.56 


67.49 


- 3.00 


67.38 


36.26 


67.34 


47.63 


67.31 


61.86 


67.30 


80.55 


67.32 


50.91 


67.30 


18. 35 


67.23 


34.54 


67.18 


53.51 


67.16 


24.04 


67.10 


26.40 


67.05 


73.72 


67.06 


77.34 


67.07 


77.93 


67.09 


124.45 


67.17 


26.01 


67.11 


52.99 


67.09 


108.13 


67.15 


91.02 


67.18 


38.62 


67.14 


154.07 


67.26 


35.74 


67.22 


100.89 


67.26 


9.30 


67.18 


157.99 


67.31 


41.82 


67.27 



725 

726 

727 

728 

729 

730 

731 

732 

733 

734 

735 

736 

737 

738 

739 

740 

741 

742 

743 

744 

745 

746 

747 

748 

749 

750 

751 

752 

753 

754 

755 

756 

757 

758 

759 

760 

761 

762 

763 

764 

765 

766 

767 

768 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 

780 

781 

782 

783 

784 

785 

786 
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3 65.96 .88 

0 0.00 0.00 

2 76.36 .66 

1 8.99 .74 

4 106.85 .95 

2 93.35 1.12 

3 106.67 .94 

3 49.35 .89 

2 73.45 .83 

2 58.83 .95 

5 160.54 .71 

3 70.26 .62 

2 145.34 .98 

2 90.10 .97 

3 175.07 .75 

3 48.96 .75 

2 39.34 1.02 

4 82.56 l . o2 

5 147.74 .95 

3 88.55 .67 

2 63.70 .77 

2 17.27 .90 

3 49.12 .93 

0 0.00 0.00 

4 108.80 .72 

6 107.57 .67 

3 35.90 .73 

2 20.87 .63 

5 133.84 1.08 

5 186.90 .86 

3 64.72 .67 

5 99.22 1.12 

4 127.06 1.06 

5 213.03 .93 

7 224.39 .83 

2 100.55 .89 

5 150.55 1.05 

1 55.21 .68 

2 161.29 .85 

4 165.35 .75 

3 66.71 .78 

1 11.03 .95 

3 29.95 .83 

3 50.62 .99 

1 14.84 1.05 

98.57 .95 

241.27 .70 

102.70 1.06 

179.70 .66 

25.59 .98 

139.35 .77 

51.32 1.05 

126.96 .93 

98.31 1.08 

5.87 .97 

49.56 .74 

125.36 .82 

97.79 .74 

26.65 .77 

53.65 .76 

48.13 .72 

50.09 .92 



55.81 


67.26 


- 3.00 


67.16 


47.83 


67.13 


3.81 


67.05 


99.31 


67.09 


102.01 


67.14 


98.06 


67.18 


41.37 


67.15 


58.36 


67.13 


53.34 


67.12 


111.90 


67.18 


41.19 


67.14 


140.01 


67.24 


85.12 


67.26 


128.82 


67.35 


33.94 


67.30 


37.36 


67.26 


81.44 


67.28 


137.67 


67.38 


57.21 


67.36 


46.53 


67.33 


12.78 


67.26 


43.22 


67.23 


- 3.00 


67.13 


75.85 


67.15 


70.47 


67.15 


23.64 


67.09 


10.44 


67.02 


141.96 


67.12 


158.81 


67.24 


40.67 


67.20 


108.39 


67.26 


132.54 


67.34 


195.71 


67.51 


184.26 


67.67 


87.13 


67.69 


156.19 


67.81 


34.80 


67.77 


134.61 


67.85 


121.49 


67.92 


49.30 


67.90 


7.65 


67.82 


22.37 


67.76 


47.72 


67.74 


12.71 


67.66 


91 .34 


67.69 


166.65 


67.82 


105.95 


67.87 


116.54 


67.94 


22.49 


67.88 


105.24 


67.92 


51.30 


67.90 


115.05 


67.96 


103.50 


68.01 


2.84 


67.93 


34.43 


67.88 


99.70 


67.92 


70.35 


67.93 


17.77 


67.86 


38.26 


67 ,. 83 


32.15 


67.78 


43.57 


67.75 



787 

788 

789 

790 

791 

792 

793 

794 

795 

796 

797 

798 

799 

800 

801 

802 

803 

804 

805 

806 

807 

808 

809 

810 

811 

812 

813 

814 

815 

816 

817 

818 

819 

820 

821 

822 

823 

824 

825 

826 

827 

828 

829 

830 

831 

832 

833 

834 

835 

836 

837 

838 

839 

840 

841 

842 

843 

844 

845 

846 

847 

848 
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268.35 


.69 


183.17 


67.90 


52.16 


.86 


42.73 


67.86 


227.35 


.58 


129.15 


67.94 


123.37 


.70 


83.23 


67.96 


220.95 


.99 


216.57 


68.15 


77.59 


.88 


65.76 


68 . 15 


0.00 


0.00 


- 3.00 


68.06 


62.05 


.85 


49.93 


68.03 


149.27 


.67 


97.15 


68.07 


64.55 


.72 


44 . 05 


68.04 


196.83 


.79 


154.07 


68.15 


39.87 


1.02 


30.02 


68.11 


121.40 


.65 


76.52 


68.12 


81.34 


.91 


71.68 


68.12 


0.00 


0.00 


- 3.00 


68.04 


64.46 


.78 


48.03 


68.01 


35.98 


1.05 


35.22 


67.97 


32.20 


.72 


20.47 


67.91 


59.10 


.78 


43.69 


67.88 


107.40 


.58 


59.73 


67.87 


100.10 


.83 


80.88 


67.89 


42.48 


.74 


28.83 


67.84 


22.94 


.72 


13.76 


67.77 


143.80 


.98 


139.14 


67.86 


39.39 


.78 


20.00 


67.81 


121.03 


.84 


99.36 


67.85 


7.69 


.89 


4.03 


67.77 


102.77 


.86 


86.06 


67.79 


150.14 


.89 


130.44 


67.87 


53.61 


.53 


25.83 


67.82 


101.02 


.91 


89.00 


67.84 


129.17 


1.4 


131.75 


67.92 


41.61 


.84 


32.13 


67.88 


103.07 


.84 


84.30 


67.90 


112.40 


.67 


73.00 


67.91 


135.28 


.86 


114.57 


67.96 


148.40 


.84 


123.12 


68.03 


45.92 


1.01 


43.66 


68.00 


147.20 


.84 


121.28 


68.06 


160.11 


.96 


152.26 


68.17 


310.18 


.93 


286.23 


68.43 


82.17 


.98 


77.82 


68.44 


115.43 


.79 


89.25 


68.47 


10.69 


.95 


7.29 


68.39 


143.61 


1.09 


154.52 


68.50 


121.62 


.74 


87.14 


68.52 


141.98 


.82 


113.49 


68.57 


153.20 


.76 


113.75 


68.63 


0.00 


0.00 


- 3.00 


68.54 


214.81 


.91 


193.06 


68.69 


16.05 


.86 


11.00 


68.62 


80.43 


.79 


60.80 


68.61 


26.99 


.59 


12.99 


68.55 


35.92 


.81 


26.60 


68.50 


119.20 


.90 


105.59 


68.54 


17.49 


.56 


7.15 


68.47 


99.31 


1.07 


103.34 


68.51 


15.99 


.97 


12.70 


68.44 


133.08 


.85 


111.29 


68.49 


80.97 


.87 


68.18 


68.49 


16.01 


.99 


13.22 


68.43 


75.06 


.65 


46 . 46 ' 


68.40 



849 

850 

851 

852 

853 

854 

855 

856 

857 

858 

859 

860 

861 

862 

863 

864 

865 

866 

867 

868 

869 

870 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 

881 

882 

883 

884 

885 

886 

887 

888 

889 

890 

891 

892 

893 

894 

895 

896 

897 

898 

899 

900 

901 

902 

903 

904 

905 

906 

907 

908 

909 

910 
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76.49 


.88 


65.20 


68.40 


131.13 


1 .08 


139.69 


68.48 


55.97 


1.02 


54.67 


68.47 


175.05 


.85 


145.97 


68.56 


148.64 


.64 


93.43 


68.59 


155.51 


.75 


114.18 


68.64 


99.10 


1.02 


98.70 


68.67 


61.57 


.68 


39.06 


68.64 


45.98 


.79 


33.65 


68.60 


42.86 


.77 


30.35 


68.55 


57.38 


.73 


39.48 


68.52 


278.61 


.80 


219.98 


68.70 


7.82 


.71 


2.71 


68.62 


21.67 


.72 


12.95 


68.55 


89.49 


.97 


84.41 


68.57 


132.26 


.90 


115.94 


68.63 


16.65 


.88 


11.73 


68.56 


75.61 


.65 


46.41 


68.54 


30.31 


.63 


16.38 


68.48 


59.10 


.90 


50.14 


68.46 


51.81 


.80 


38.95 


68.42 


78.80 


.82 


62.07 


68.41 


35.14 


.85 


27.39 


68.37 


42.67 


.89 


35.67 


68.33 


92.72 


.69 


61.14 


68.32 


25.99 


1.09 


25.58 


68.27 


7.83 


.55 


1.49 


68.20 


98.31 


.79 


74.79 


68.20 


60.85 


.72 


41.18 


68.17 


82.45 


.91 


72.05 


68.18 


46.22 


.64 


27.37 


68 . 13 


105.05 


1.04 


107.74 


68.18 


72.59 


.61 


41.54 


68.15 


131.70 


.98 


126.04 


68.21 


42.01 


1.03 


40.67 


68.18 


55.18 


.95 


49.75 


68.16 


49.31 


.62 


27.65 


68.11 


33.63 


1.07 


33.22 


68.07 


123.45 


.84 


101.33 


68.11 


42.31 


1.12 


44.73 


68.09 


133.15 


.91 


119.60 


68.14 


118.14 


.90 


104.31 


68.18 


117.93 


.99 


113.85 


68.23 


175.85 


.88 


153.03 


68.33 


89.84 


.78 


67.39 


68.33 


14.98 


1.04 


12.70 


68.27 


60.13 


.79 


45.10 


68.24 


76.65 


.60 


43.22 


68.21 


36.35 


.92 


30.89 


68.17 


18.17 


1.11 


17.27 


68.11 


60.94 


.94 


54.87 


68.10 


69.61 


.96 


63.87 


68.10 


165.45 


.61 


98.87 


68.13 


44.09 


.79 


32.31 


68.09 


121.40 


.81 


95.86 


68.12 


32.84 


.99 


29.82 


68.08 


88.36 


.91 


77.91 


68.09 


6.56 


.65 


1.42 


68.02 


299.06 


.92 


274.07 


68.24 


67.80 


.90 


58.48 


68.23 


98.84 


.71 


67.79 


68.23 


64.58 


.96 


59.49 


68.22 



911 

912 

913 

914 

915 

916 

917 

918 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 

931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 

943 

944 

945 

946 

947 

948 

949 

950 

951 

952 

953 

954 

955 

956 

957 

958 

959 

960 

961 

962 

963 

964 

965 

966 

967 

968 

969 

970 

971 

972 
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2 64,88 1.01 

75.28 1.11 

120.56 .95 

103.66 .79 

59.12 .90 

106.89 .95 

117.69 .98 

91.22 .96 

120.53 .68 

71.18 .69 

63.55 .97 

41.41 .79 

61.95 .74 

1 17.27 .64 

2 28.94 1.04 

3 57.08 1.03 

4 160.68 .93 

2 44.95 .86 

4 57.39 .67 

6 285.64 .76 

2 18.85 .76 

5 131.95 1 -01 

1 24.78 1.05 

4 149.44 .82 

2 25.00 .83 

2 108.64 .81 

6 127.18 .60 

1 34.01 1.01 

3 68.04 .76 

4 114.78 .83 

3 93.84 .74 

1 12.16 .90 

6 235.57 .91 

5 141.75 .90 

4 86.68 .91 

1 95.48 .83 

1 12.50 1.00 

3 94.18 .82 

3 131.53 .81 

2 46.95 .94 

0 0.00 0.00 

2 41.66 1.00 

5 107.98 .69 

3 112.16 .81 

2 58.17 1.00 

1 18.07 .73 

4 226.90 .62 

3 71.78 .93 

2 18.57 .81 

i 2 120.67 1.09 

2 74.21 .80 

2 9.77 .94 

3 102.85 .81 

2 48.64 .67 

4 134.03 .88 

4 141.04 1.07 

2 58.36 .94 

3 43.71 .90 

4 245.90 .76 

4 118.89 .92 

3 81.39 .95 

3 99.31 -67 



62.52 


68.22 


81.36 


68.23 


112.32 


68.28 


79.48 


68.29 


50.52 


68.27 


98.41 


68.30 


112.95 


68.35 


84.89 


68.37 


78.87 


68.38 


46.95 


68.36 


58.82 


68.35 


30.41 


68.31 


43.02 


68.28 


8.23 


68.21 


27.33 


68.17 


56.27 


68.16 


146.55 


68.24 


35.79 


68.21 


36.20 


68.17 


213.59 


68.33 


11.59 


68.27 


131.12 


68.34 


23.11 


68.29 


119.74 


68.34 


18.00 


68.29 


85.45 


68.31 


73.65 


68.31 


31.54 


68.27 


49.10 


68.25 


92.66 


68.28 


66.81 


68.28 


8.12 


68.21 


212.77 


68.37 


126.02 


68.43 


76.62 


68.44 


76.82 


68.45 


9.69 


68.38 


74.29 


68.39 


104.49 


68.43 


41.47 


68.40 


- 3.00 


68.32 


39.17 


68.29 


72.56 


68.30 


88.83 


68.32 


55.31 


68.31 


10.31 


68.25 


139.19 


68.32 


64.04 


68.31 


12.40 


68.26 


128.42 


68.32 


56.32 


68.31 


6.46 


68.24 


80.43 


68.26 


29.81 


68.22 


115.00 


68.26 


148.36 


68.35 


52.05 


68.33 


36.96 


68.30 


184.27 


68.42 


106.62 


68.46 


74.41 


68.46 


63.99 


68.46 



973 


3 


75.64 


1.01 


74.07 


68.46 


974 


5 


68.62 


.97 


64.64 


68.46 


975 


2 


49.27 


.83 


37.95 


68.43 


976 


4 


120.25 


.97 


114.25 


68.48 


977 


4 


113.96 


.71 


78.62 


68.49 


978 


1 


9.93 


.97 


6.76 


68.42 


979 


1 


23.98 


.83 


17.12 


68.37 


980 


0 


0.00 


0.00 


-3.00 


68.30 


981 


3 


71.60 


• 80 


54.54 


68.28 


982 


3 


36.94 


.91 


30.93 


68.25 


983 


2 


48.01 


.67 


29.35 


68.21 


984 


2 


59.59 


.74 


41.24 


68.18 


985 


2 


91.15 


.59 


50.75 


68.16 


986 


1 


7.20 


1.10 


5.07 


68.10 


987 


3 


89.59 


.80 


69.07 


68.10 


988 


1 


33.57 


.63 


18.46 


68.05 


989 


2 


50.60 


1.07 


51.47 


68.03 


990 


2 


96.38 


.74 


68.50 


68.03 


991 


4 


47.68 


.82 


36.88 


68.00 


992 


3 


91.39 


1.01 


89.55 


68.02 


993 


3 


30.18 


.97 


26.77 


67.98 


994 


7 


170.25 


.82 


138.48 


68.05 


995 


3 


54.44 


.77 


39.36 


68.02 


996 


8 


237.43 


.73 


172.56 


68.13 


997 


1 


13.36 


.99 


10.38 


68.07 


998 


3 


61.66 


.58 


33.07 


68.03 


999 


3 


112.05 


.90 


98.59 


68.06 


1000 


3 


43.11 


.63 


24.54 


68.02 


THE BEST 


ESTIMATE OF 


PROFIT BY 


THE MONTE 


CARLO SIMULATION 


METHOD 


AT 1000 


iterations 


IS 68.02 


MILLION 1 


DOLLARS 





THE VARIANCE IS 2879.367 

THE STANDARD DEVIATION IS 53.660 

THE COEFFICIENT OF VARIATION IS .7889 



PROFIT DISTRIBUTION OF AN EXPLORATORY DRILLING PROGRAM 



NO. OF WELLS 
DRILLED 
20 



PROB OF DRY HOLE 
SUCCESS COST QMEDIAN Q04 
.150 .150 20.00 46.00 



QMAX 

140.00 



QMIN 

1.00 



THE TRIANGULAR DISTRIBUTION FOR DETERMINATION OF THE 
PROFIT PER BARREL HAS THE FOLLOWING EXTREMES IN CENTS - 
A = 50.00 B = 88.00 C = 115.00 



UPPER LIMIT OF PROBABILITY CUMMULATIVE 

PROFIT INTERVAL OF interval probability 



-10 


*0000 


*0000 


0 


.0420 


.0420 


10 


.0530 


.0950 


20 


.0740 


.1690 


30 


.0910 


.2600 


40 


• 1060 


.3660 


50 


.0870 


.4530 


60 


.0780 


.5310 


70 


.0630 


.5940 


80 


.0630 


.6570 


90 


.0650 


.7220 


100 


.0450 


.7670 


110 


.0490 


.8160 


120 


.0340 


.8500 


130 


.0280 


.8780 


140 


.0280 


.9060 


150 


.0130 


.9190 


160 


.0220 


.9410 


170 


.0090 


.9500 


180 


.0050 


.9550 


190 


.0110 


.9660 


200 


.0070 


.9730 


210 


.0050 


.9780 


220 


.0050 


.9830 


230 


.0030 


.9860 


240 


.0030 


.9890 


250 


*0000 


.9890 


260 


.0040 


.9930 


270 


.0020 


.9950 


280 


.0030 


.9980 


290 


.0010 


.9990 


300 


*0000 


.9990 


350 


*0000 


.9990 


400 


.0010 


1.0000 


450 


*0000 


1.0000 


500 


*0000 


1.0000 


550 


*0000 


1.0000 
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